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EXECUTIVE  SUMMARY 


A  self-consistent  materials  chemistry  of  flaming  combustion  is  derived  from  solid-state  kinetics  and 
thermodynamics  using  a  realistic  physical  model  of  polymer  burning.  The  rates  of  solid-state 
(pyrolysis)  and  gas  phase  (combustion)  reactions  are  assumed  to  be  rapid  in  comparison  to  the  rate 
of  heat  transfer  at  the  burning  surface.  Coupling  of  thermal  diffusion  and  chemical  kinetics  occurs 
in  the  surface  pyrolysis  zone  where  the  rate  of  temperature  rise  for  steady  burning  becomes  the 
characteristic  heating  rate  for  thermal  decomposition. 

Detailed  thermal  degradation  chemistry  is  foregone  in  favor  of  a  transient  mass  balance  on  the 
polymer,  fuel  gases,  and  solid  char  in  the  anaerobic  pyrolysis  zone.  Closed-form,  time-independent 
solutions  for  the  scalar  mass  loss  rate  and  char  yield  are  obtained  from  the  degradation  kinetics 
which,  in  combination  with  the  solid-state  thermal  transport  and  thermodynamic  properties 
calculated  from  the  polymer  chemical  structure,  provide  the  scaling  relationship  between  material 
properties  and  steady-burning  rate.  The  critical  condition  for  a  nonzero  burning  rate  in  the  absence 
of  an  external  heat  flux  (flammability)  can  then  be  expressed  in  terms  of  material  properties. 

Thermochemical  predictions  of  ignitability  and  flaming  heat  release  rate  for  a  variety  of  charring 
and  noncharring  polymers  are  in  general  agreement  with  experimental  data  using  the  present 
approach. 
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1.  INTRODUCTION. 
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Synthetic  polymers  have  improved  our  quality  of  life  since  their  commercial  introduction  over  50 
years  ago.  Synthetic  polymers  can  be  low-cost  functional  replacements  for  metals  or  natural 
polymers  (e.g.,  wool,  silk,  leather,  and  wood)  with  improved  physical,  chemical,  and  optical 
properties  or  can  be  enabling  materials  for  high-technology  applications  such  as  electronics,  optics, 
medical  devices,  and  aerospace.  Commodity  applications  of  synthetic  polymers  include  home 
furnishings,  appliances,  construction,  recreation,  and  consumer  packaging;  but  it  is  the  performance- 
critical,  high-technology  applications  that  have  driven  the  development  of  thermally  and  chemically 
stable  polymers  with  intrinsic  fire  resistance  [1,2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12]. 

The  use  of  polymers  in  public  buildings  and  mass  transportation  where  fire  is  a  significant  threat  to 
human  life  because  of  the  potential  for  rapid  fire  spread  calls  for  a  balance  between  the  functionality 
of  polymers  and  their  risk  as  a  fire  hazard.  Flames  from  fires  destroy  property,  but  it  is  the 
combustion  products  in  the  smoke  which,  directly  or  indirectly,  lead  to  loss  of  life  [13,  14,  15]. 
Conventional  methods  of  reducing  ignitability  in  small  flame  tests  by  the  addition  of  flame- 
retarding  chemicals  usually  increases  the  density  and  toxicity  of  smoke  and  fumes  [13,  16]. 
Although  flame-retardant  materials  may  bum  less  vigorously  under  real  fire  conditions  [16],  the 
lowered  bum  rate  is  usually  accompanied  by  higher  smoke  levels  resulting  from  incomplete  gas 
phase  combustion.  Incomplete  combustion  increases  the  carbon  monoxide/carbon  dioxide  ratio 
and  the  toxic  potency  of  the  smoke  [14, 16],  Thus,  the  most  effective  way  to  improve  the  fire  safety 
of  components  is  to  use  intrinsically  fire-resistant  polymeric  materials  [17,  18,  19]  with  reduced 
burning  rate  and  high  gas  phase  combustion  efficiency. 

The  above  and  subsequent  references  to  polymer  burning  are  understood  to  mean  the  gas  phase 
combustion  of  the  mixture  of  oxygen  and  polymer  thermal  degradation  products  in  a  fire,  since 
technically  polymers  do  not  bum  in  the  condensed  (solid  or  liquid)  state.  Polymers  do  not  bum  in 
the  condensed  state  because  of  the  low  solubility  and  diffusivity  of  the  oxidizer  (diatomic  oxygen) 
and  the  low  oxidation  rate  at  the  decomposition  temperature.  In  fact,  thermal  degradation  of  the 
surface  layer  of  polymer  in  the  presence  of  a  heat  source  is  thought  to  occur  in  a  reducing,  rather 
than  an  oxidizing,  environment.  Low  molecular  weight  volatile  organic  compounds  are  produced 
which  mix  with  atmospheric  oxygen  above  the  polymer  surface  to  form  a  flammable  mixture  which 
when  ignited  combusts,  producing  a  luminous  flame.  Of  particular  interest  is  the  rate  of  fuel 
generation  at  the  polymer  surface  in  a  fire  since  this  process  governs  the  heat  release  rate  and 
smoke  evolution  during  gas  phase  combustion. 

Figure  1  illustrates  the  three  coupled  processes  required  for  flaming  combustion:  (1)  heating  of  the 
polymer,  (2)  thermal  decomposition  to  gaseous  products,  and  (3)  ignition  and  combustion  in  air. 
An  ignition  source  or  thermal  feedback  of  radiant  energy  from  die  flame  supplies  heat  to  the 
polymer  surface  which  causes  thermolytic  cleavage  of  primary  chemical  bonds  in  the  polymer 
molecules.  Evaporation  of  the  low  molecular  weight  degradation  products  and  reaction  with  air 
(oxygen)  in  the  combustion  zone  above  the  surface  releases  heat  and  produces  carbon  dioxide, 
water,  and  incomplete  combustion  products  such  as  carbon  monoxide,  unbumed  hydrocarbons,  and 
soot.  In  order  to  resist  burning,  the  fire  cycle  must  be  broken  at  one  or  more  places. 
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FIGURE  1.  SCHEMATIC  DIAGRAM  OF  POLYMER  BURNING 

Whether  a  polymer  is  fire  resistant  or  contributes  to  the  spread  of  fire  depends  on  the  intensity  of 
the  fire  or  heat  source  and  the  polymer’s  reaction  to  it  (e.g.,  fuel  generation  rate).  Heat  fluxes  range 
in  intensity  from  a  small  ignition  source  such  as  a  match  or  cigarette  20  W),  to  a  Bunsen  burner 
flame  (=  102  W),  to  a  trash  can  fire  (~  105  W),  to  a  full-sized  room  fire  (=  106  W),  and  to  a 
building  or  hydrocarbon  pool  fire  (~  107  W)  [20,  21].  Surprisingly,  the  heat  flux,  which  is  the  rate 
at  which  heat  enters  per  unit  area  of  an  adjacent  material  surface,  falls  into  the  relatively  narrow 
range  [20]  of  20-200  kW/m2  for  this  107  range  of  fire  sizes.  It  is  the  area  of  the  heated  surface, 
not  the  heat  flux  itself,  which  is  proportional  to  the  fire  size  [20]. 

Figure  2  shows  the  maximum  possible  surface  temperature  for  an  £  =  0.75  emissivity  material  as  a 
function  of  incident  heat  flux  for  a  surface  convective  heat  transfer  coefficient  h  =  15  W/m2-K 
assuming  thermal  equilibrium  [22] 

(W  =  eo[Ts4-T04]  +  fi[Ts-T0]  (1) 

Ranges  of  heat  flux  incident  on  a  material  in  proximity  to  fires  of  different  sizes  are  shown  parallel 
to  the  horizontal  (heat  flux)  axis  in  figure  2.  The  vertical  axis  is  subdivided  into  qualitative  fire 
performance  rankings  which  presume  short-term  thermal  stability  (decomposition  temperature)  in 
the  range  of  surface  temperatures  experienced  by  a  material  in  various  fire  environments.  Small 
diffusion  flames  of  burning  material  have  heat  fluxes  in  the  range  of  25-40  kW/m2  corresponding 
to  the  range  of  sample  surface  temperatures  500-600°C  in  figure  2.  Polymers  with  short-term 
thermal  stability  in  this  range  of  temperatures  are  typically  flame  resistant  since  they  will  not 
continue  to  bum  in  the  absence  of  an  external  heat  source  or  elevated  oxygen  concentration.  Small 
fires  in  proximity  to  the  polymer  surface  (e.g.,  a  burning  wastebasket)  or  remote  large  fires  generate 
surface  heat  fluxes  in  the  range  30-60  kW/m2  depending  on  their  size  and  distance  from  the 
surface  Heat  fluxes  of  30-60  kW/m2  generate  equilibrium  surface  temperatures  in  the  550-750°C 
range.  Polymers  with  short-term  thermal  stability  in  this  range  of  temperatures  may  be  categorized 
as  fire  resistant  since  they  typically  resist  piloted  ignition  or,  when  ignited,  bum  with  a  low  rate  of 
heat  release  (<  100  kW/m2).  Intrinsically  fire-resistant  polymers  are  thermally  stable  polycyclic 
and  heteroatomic  molecules  with  low  hydrogen  content.  Noncombustible  behavior  in  figure  2  and 
in  standard  tests  [23]  are  associated  with  short-term  thermal  stability  at  or  above  temperatures  of 
750°C  corresponding  to  heat  fluxes  in  excess  of  75  kW/m2  in  figure  2.  These  large  heat  fluxes 
result  from  contact  with,  or  close  proximity  to,  fully  developed  large  fires  such  as  liquid 
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hydrocarbon  pool  fires,  burning  buildings,  or  post-flashover  compartments.  Noncombustible 
materials  are  presently  inorganic  in  composition  [24] . 


FIGURE  2.  EQUILIBRIUM  SURFACE  TEMPERATURE  OF  POLYMER  VERSUS 
EXTERNAL  HEAT  FLUX  FOR  FIRES  OF  VARIOUS  SIZES 

Several  comprehensive  texts  have  been  written  on  the  gas  phase  processes  of  flaming  combustion 
[25, 26, 27, 28]  while  the  thermochemical  processes  which  occur  in  the  solid  state  to  generate  the 
combustible  gases  in  a  fire  have  received  relatively  little  attention  [29, 30].  The  remainder  of  this 
chapter  explores  the  material  science  of  flaming  combustion  by  developing  relationships  between 
the  chemical  structure  of  polymers  and  their  burning  behavior.  Recent  developments  in  the 
metrology  and  modeling  of  fire  [31,  32]  and  its  impact  on  materials  [29,  33,  34,  35]  provide  a 
physical  basis  for  polymer  ignition  and  burning  in  terms  of  measurable,  macroscopic  flammability 
parameters.  Connecting  these  macroscopic  flammability  parameters  to  the  molecular  structure  of 
the  polymer  through  the  kinetics  and  thermodynamics  of  the  fuel  generation  process  provides  a 
thermochemical  basis  for  the  solid-state  processes  of  flaming  combustion. 

2.  FIRE  BEHAVIOR  OF  COMBUSTIBLE  SOLIDS. 

The  solid-state  processes  in  the  flaming  combustion  of  polymers  which  can  be  treated  at  the 
continuum  level  are  the  subject  of  this  section.  The  continuum  treatment  disregards  the  discrete 
(molecular)  structure  of  matter  so  that  the  temperature  distribution  and  more  importantly  its 
derivatives  are  continuous  throughout  the  material.  In  addition,  the  material  is  assumed  to  have 
identical  thermal  properties  at  afi  points  (homogeneous)  and  in  all  directions  (isotropic).  The 
concept  of  a  continuous  medium  allows  us  to  define  fluxes  at  a  point,  e.g.,  a  surface  in  one¬ 
dimensional  space.  Chemical  reactions  in  the  solid  (pyrolysis)  and  flame  (combustion)  are 
assumed  to  occur  so  rapidly  that  the  burning  rate  is  determined  solely  by  the  heat  transfer  rate  [26], 
Differential  [36,  37,  38,  39]  and  integral  [40,  41]  condensed-phase  burning  models  have  been 
developed  from  the  continuum  perspective  with  coupled  heat  and  mass  transfer  for  both  charring 
[36,  37,  40]  and  noncharring  polymers  [39,  41].  All  of  these  models  must  be  solved  numerically 
for  the  transient  (time-dependent)  mass  loss  and  heat  release  rates. 
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In  the  present  treatment  we  neglect  transient  phenomena  in  thick  samples,  such  as  char  layer  growth 
and  the  progression  of  crystalline  melting,  and  consider  only  the  quasi-steady  burning  of  an 
idealized  solid  with  constant  surface  heat  flux.  This  simplified  approach  allows  us  to  develop 
algebraic  (scaling)  relationships  between  material  composition  and  the  solid-state  processes  of 
flaming  combustion  but  ignores  many  important  details  of  thick-sample  burning  behavior  which 
can  only  be  captured  through  transient  numerical  analysis. 


2.1  IGNITION. 

Figure  3  shows  the  geometry  of  the  continuum  model  of  polymer  combustion.  If  the  polymer 
thickness  b+8  in  figure  3  is  large  compared  to  the  thermal  degradation  process  (pyrolysis)  zone 
thickness  8  (i.e.,  b  »  8),  then  according  to  the  critical  surface  temperature  criteria  for  ignition, 
sustained  piloted  ignition  occurs  when  the  polymer  surface  reaches  its  thermal  decomposition 
temperature.  A  critical  surface  temperature  is  only  one  of  several  criteria  which  have  been  proposed 
for  ignition  [42],  Other  critical  values  at  ignition  include  the  average  temperature  of  the  solid, 
pyrolyzate  mass  flux,  char  depth,  rate  of  increase  of  local  gas  temperature,  reaction  rate  in  the  gas 
boundary  layer,  and  a  gas  temperature  gradient  reversal.  As  will  be  shown  in  the  following  sections, 
the  decomposition  temperature  of  a  polymer  is  a  kinetic  parameter  whose  exact  value  depends  on 
the  heating  rate  (heat  flux). 


Heat  transfer  in  the  solid  polymer  is  described  by  the  one-dimensional  energy  equation  for 
unsteady  heat  conduction  with  no  internal  heat  generation 


3T  ... 
W  ~ 


(2) 


where  T  is  the  temperature  at  location  x  in  the  solid  polymer  and  a  =  K/pc  is  the  polymer  thermal 
diffusivity  in  terms  of  its  thermal  conductivity,  K,  density,  p,  and  heat  capacity,  c;  and  v  is  the 
regression  velocity  of  the  burning  surface.  During  the  preheat  phase  prior  to  ignition,  there  is  no 
surface  regression,  so  v  =  0  and  equation  2  reduces  to 


82T  1  8T 

9x2  a  dt 


(3) 


The  solution  [43]  of  equation  3  for  the  ignition  time  tign  of  a  thermally  thick  sample  with  a  constant 
net  heat  flux  qnet  at  the  surface  x  =  0  is 


vign 


n 

4  KPC 


Qnet 


where  %  is  the  (piloted)  ignition  temperature  which  is  approximately  equal  to  the  peak  mass  loss 
temperature  T  for  transient  heating  [44],  and  T0  is  the  ambient  initial  temperature.  If  8  is  greater 
than  the  sampfe  thickness  the  sample  is  considered  thermally  thin,  and  ignition  occurs  at  time 


tig„  =  Pbc 
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Equations  4  and  5  state  that  the  time  to  ignition  is  determined  by  the  ignition  (decomposition) 
temperature;  the  material’s  thermal  and  transport  properties  K,  p,  and  c;  and  the  net  heat  flux  to  the 
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surface.  Equations  4  and  5  are  accurate  when  the  incident  heat  flux  is  high  compared  to  heat  losses 
by  surface  convection  and  reradiation. 


FIGURE  3.  GEOMETRY  OF  POLYMER  COMBUSTION  ANALYSIS 
2.2  STEADY  BURNING. 

Once  sustainable  ignition  has  occurred,  steady,  one-dimensional  burning  of  the  polymer  is 
assumed.  Steady  burning  at  a  constant  surface  heat  flux  is  treated  as  a  stationary  state  by  choosing 
a  coordinate  system  which  is  fixed  to  the  surface  and  moving  at  the  recession  velocity  v.  If  there  is 
no  internal  heat  generation  or  absorption,  the  one-dimensional  heat  conduction  equation 
(equation  2)  applies.  Since  semicrystalline  polymers  absorb  the  heat  of  fusion  during  melting  at 
temperatures  below  the  decomposition  temperature,  equation  2  is  only  approximate  for  these 
materials.  Under  steady-state  conditions,  dT(x)/dt  =  0  and  the  constant  thickness  pyrolysis  zone 
moves  through  the  infinitely  thick  solid  with  a  fixed  temperature  gradient  so  that  equation  2 
becomes 


d2T  .  y  dT 
dx2  a  dx 


0 


(6) 


for  steady  burning  of  a  material  with  a  constant  thermal  diffusivity  a  =  K/pc.  The  general  solution 
of  equation  6  is 

T(x)  =  c,  +  c2exp[-vx/a]  (7) 

Two  boundary  conditions  are  needed  to  evaluate  the  constants  of  integration  ci  and  C2  in 
equation  7.  Conservation  of  energy  at  the  pyrolysis  front  x  =  0  gives 


flfT(x) 
dx  x=0 


=  -q„et  +  pvAhv 


(8) 


from  which  C2  =  (qnet  cc/kv)  -  (Ah/c)  with  Ahv  the  latent  heat  of  vaporization  of  the  pyrolysis 
products  and  qnet  the  net  heat  flux  at  the  surface  (x  =  0)  expressed  as 
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(9) 


9  net  =  9 ext  +  9 flame  “  EO(^  ^o)  ^(Ts  T0) 


—  Qext  ^  Aflame  4cr 

Equation  9  defines  the  net  heat  flux  into  the  surface  qnet  as  the  difference  between  the  heat  flux 
entering  the  surface  from  an  external  radiant  energy  source,  qext ,  and/or  surface  flame,  9flame ,  and 
the  critical  heat  flux  for  ignition,  qcr .  The  critical  heat  flux  for  ignition  is  equal  to  the  heat  removed 
by  reradiation,  ec(Ts4  -  T04),  and  convection,  h(Ts  -  T0),  at  the  ignition  temperature  Ts  =  Tign  =  Tp 
where  e,  a,  and  h  are  the  surface  emissivity,  Stefan-Boltzmann  constant,  and  average  surface 
convective  heat  transfer  coefficient,  respectively.  In  the  absence  of  an  external  radiant  energy  source 
(e.g.,  during  a  Bunsen  burner-type  flammability  test),  qnet  =  0  when  qcr  =  9flame  • 

On  the  rear  face  of  the  infinite  slab  (x  =  °°),  we  specify  dT/dx  =  0  or  equivalently  T(°°)  =  T0  =  c, 
where  T0  is  the  ambient  temperature  in  equation  8.  The  final  temperature  distribution  during 
steady-state  burning  of  an  infinitely  thick  material  is 


T(x)  =  T0  + 


9net  1 

- exp 

pc  v 


(10) 


The  steady  recession  velocity  of  the  surface  x  =  0  at  temperature  T(0)  —  Tp  from  equation  10  is 


_  1  _ 9  net  _  _  1  9  net 

p(c(Tp-T0)  +  Ahv)  "  P  hg 

where  the  total  heat  of  gasification  hg  per  unit  original  mass  of  polymer  is  [45] 

hg  =  c(Tp  -  T0)+  Ahv.  (12) 

Equations  10  and  1 1  allow  the  steady-state  temperature  distribution  in  the  burning  solid  polymer  to 
be  expressed  as 
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which  is  in  qualitative  agreement  with  experimental  data  for  the  temperature  gradient  in  steadily 
burning  liquid  pools  [32]  if  Tp  is  taken  as  the  boiling  temperature  of  the  liquid  fuel. 

Conservation  of  mass  for  the  control  volume  in  which  the  virgin  polymer  of  density  p  pyrolyzes  to 
an  inert  fraction  or  char  residue  p  gives  (personal  communication,  J.  Quintiere,  University  of 
Maryland,  1999). 
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(13) 


where  mg  is  the  mass  loss  rate  of  pyrolysis  gases  per  unit  surface  area.  Defining  a  heat  of 
gasification  per  unit  mass  of  volatiles 
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and  combining  equations  1 1  and  13 
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(14) 


The  heat  of  gasification  per  unit  mass  of  solid  polymer  Ik  can  be  determined  from  the  reciprocal 
slope  of  a  plot  of  areal  mass  loss  rate  versus  external  heat  flux  if  the  char  yield  is  measured  after  the 
test,  since  from  equations  9  and  13 


The  reciprocal  slope  of  m ^versus  qext  equals  Lg  only  if  qflame  does  not  change  with  qextor  if 
9  flame  «  q  ext  •  Multiplying  equation  15  by  the  net  chemical  heat  of  complete  combustion  of  the 
volatile  polymer  decomposition  products  h°  and  the  gas  phase  combustion  efficiency  %  (see 
section  4.3)  gives  the  usual  result  for  the  average  heat  release  rate  of  a  burning  specimen  [32, 46]. 

h°c  .  h°  . 

qc  =  %hcmg  =  X  (1  -  p)  p  q„e,  =  %jj-q„et  (16) 

Given  that  %  lies  in  the  relatively  narrow  range  [46],  X  =  0.5-0.9,  it  is  the  combustibility  ratio  [47]  or 
heat  release  parameter  [46]  h°  /Lg  =  (l-p)/hg  which  is  the  dominant  material  burning  parameter 
since  it  can  vary  by  more  than  an  order  of  magnitude  for  polymer  solids  [46].  Since  h^  is  relatively 
constant  for  a  wide  range  of  materials  (sections  4.2  and  4.3),  it  is  the  variation  in  char  yield 
p  =  0  -  0.8  and  h°  =3-45  kJ/g  which  accounts  for  most  of  the  difference  in  the  heat  release 
parameter  between  materials. 

Figure  4  shows  idealized  heat  release  rate  histories  for  steady  burning  of  thermally  thick  polymers 
according  to  the  equations  developed  in  this  section.  The  upper  curves  represent  steady  burning  of 
a  noncharring  (p  =  0)  and  charring  (p  =  0.5)  polymer,  respectively,  with  comparable  hg.  The  time- 
independent  heat  release  rate  curve  for  the  charring  material  (dashed  line)  is  hypothetical  since  the 
steady-burning  model  assumes  a  constant  surface  recession  rate,  i.e.,  no  accumulation  of  char  at  the 
surface.  The  lower  curve  (solid  line)  is  a  realistic  heat  release  history  for  a  charring  polymer 
showing  transient  effects  such  as  a  peak  in  heat  release  rate  soon  after  ignition  followed  by  a 
depression  in  the  heat  release  rate  as  the  char  layer  increases  in  thickness.  The  growing  char  layer 
insulates  the  underlying  polymer  from  the  surface  heat  flux  and  acts  as  a  diffusion  barrier  to  the 
volatile  fuel.  Charring  polymers  can  be  linear  (thermoplastic)  or  crosslinked  (thermoset)  polymers 
having  amorphous  or  semiciystalline  morphologies.  The  area  under  the  heat  release  rate  curves  per 
unit  mass  of  polymer  consumed  is  the  effective  heat  of  flaming  combustion.  The  effective  heat  of 
combustion  is  determined  primarily  by  the  combustion  chemistry  in  the  flame.  Combustion 
efficiency  decreases  when  halogens  are  present,  when  soot/smoke  is  produced  in  large  yield,  or 
when  there  is  insufficient  oxygen  for  complete  combustion.  Flaming  combustion  efficiency  is 
relatively  independent  of  the  charring  tendency  of  a  polymer. 
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FIGURE  4.  IDEALIZED  HEAT  RELEASE  RATE  CURVES  FOR  STEADY  BURNING  OF 
NONCHARRING  (TOP)  AND  CHARRING  (MIDDLE)  POLYMERS  (The  lower  curve  is 

the  actual  behavior  of  charring  polymers.) 

In  the  subsequent  development  of  the  fuel  generation  rate  using  thermal  degradation  kinetics,  it  will 
be  useful  to  know  the  rate  of  temperature  rise  of  the  polymer  surface  in  the  pyrolysis  zone  as  it 
moves  through  the  solid  at  constant  velocity,  v.  From  equations  10  and  1 1  the  effective  heating  rate 
at  the  surface  is 

dT  _  d7 
dt  x=0~  V  dx 

Typically,  Ah/c(T  -T0)  ~  0.4  ±  0.1  (see  section  4.2)  so  the  heating  rate  at  the  surface  is 
approximately 

dT 

dt  x=0 
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•  2 
Qnet 


•  2 
Qnet 


2  kPc(Tp  -  Tn)  Kphs 


(17b) 


•  2 
Qnet 


1  + 


Ah„ 


KPc  (Tp  -  T0)  (  c(Tp-T0) 


(17a) 


According  to  equation  17b  the  rate  of  surface  temperature  rise  of  a  polymer  with  Tp  «  500°C 
(723  K)  and  a  typical  Kpc  =  5  x  10s  W-s-m^-K'2  experiencing  a  50  kW/m2  net  surface  heat  flux  is 

dT/dt  =  5  K/s. 

Figure  5  is  a  plot  of  net  heat  flux  versus  surface  temperature  for  polymers  at  an  external  heat  flux 
qext  =  50  kW/m2  calculated  from  equations  9  and  16.  Plotted  in  figure  5  is  the  calculated  net  heat 
flux  versus  surface  temperature  curves  for  a  horizontal  plate  with  an  average  surface  free  convective 
heat  transfer  coefficient,  h  =  15  W/m2-K,  for  polymer  emissivities  e  =  0.75  and  1.0.  Calculations 
such  as  those  plotted  in  figure  5  can  be  used  to  determine  the  ignitability  and  heat  release  rate  of  a 
polymer  as  follows.  The  minimum  (net)  heat  flux  for  ignition  can  be  calculated  using  the  criteria 

that  ignition  occurs  when  the  mass  flux  exceeds  the  critical  value  mg  ~  3-5  g-m_2-s-1  [42].  For 
typical  values  %  =  0.8  and  hc#  =  30  kJ-g"1,  the  critical  mass  flux  corresponds  to  a  critical  heat 
release  qc,CT  =  %  hc°  mg  »  100  kW/m2  (46). 
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Setting  qc  =100  kW/m2  in  equation  16,  the  critical  net  heat  flux  for  ignition  of  a  nonchaning 
polymer  with  typical  hg  =  Lg  =  3  kJ/g  is 


100kW/m2 

Xh°^g 


100kW/m2 
(0.8)(30  kJ/g)/(3  kJ/g) 


»  13  kW/m2 


which  is  in  the  range  qnet,cr  =10-15  kW/m2  measured  for  nonhalogen  commodity  polymers  [46]. 

Figure  5  shows  that  polymers  which  can  sustain  surface  temperatures  in  the  range  600-650°C 
without  thermally  degrading  to  gaseous  fuel  are  able  to  reradiate  and  convect  away  enough  heat  at 
qext  =  50  kW/m2  to  ensure  that  the  absorbed  (net)  heat  flux  is  below  the  critical  value  for  ignition 
indicated  by  the  dashed  line.  This  result  is  in  agreement  with  experimental  data  for  the  thermally 
stable  benzobisoxazole  polymer  [48]  which  decomposes  at  temperatures  above  650°C  and  is  thus 
able  to  resist  ignition  at  qext  =  50  kW/m2.  Figure  5  also  shows  that  the  emissivity  of  the  polymer 
surface  has  a  significant  effect  on  the  net  heat  flux  into  the  solid.  Low-emissivity  (heat  reflective) 
coatings  are  used  commercially  for  ignition-resistant  firefighting  apparel,  paints,  and  aircraft 
evacuation  slides. 

Improved  ignition  resistance  can  also  be  achieved  by  reducing  the  magnitude  of  the  quantity  %  h°/Lg 

so  that  the  critical  heat  flux  for  ignition  increases.  Polymer  formulators  use  this  strategy  to  pass 
small-scale  flame  resistance  tests  using  chemical  flame-retardant  additives  [49].  Halogenated 
compounds  are  routinely  added  to  polymers  to  interrupt  the  gas  phase  combustion  process  and 
reduce  %.  Metal  hydrates  which  decompose  endothermically  to  inert  gases  (e.g.,  H20)  when  heated 
are  added  to  dilute  the  combustibles  in  the  solid  and  in  the  flame  [50].  Inert  compounds  and  char 
promoters  increase  Lg  by  allowing  the  polymer  to  absorb  heat  without  liberating  fuel.  Some  flame- 
retardant  additives  work  by  a  combination  of  mechanisms  (i.e.,  synergistically),  but  the  overall  effect 
is  to  elevate  the  critical  heat  flux  for  ignition  in  the  flame  test  until  the  flame  retardant  is  depleted 
from  the  polymer. 


FIGURE  5.  NET  HEAT  FLUX  AT  qext  =  50  kW/m2  VERSUS  POLYMER  SURFACE 
TEMPERATURE  FOR  EMISSIVITIES  e  =  0.75  and  1.0  (Dashed  line  is  the  critical  heat 

flux  for  %h°  /Lg  =  8.) 


9 


The  continuum  analyses  demonstrated  the  critical  parameters  for  fire-resistant/low  heat  release  rate 
polymers:  high  decomposition  temperature  T_,  low  fuel  value  of  degradation  products  hc°,  a  large 
heat  of  gasification  h ,  and  high  char  yield  |i.  While  low  combustion  efficiency  x  is  a  viable  route 
to  ignition-resistant  materials,  the  increased  hazard  associated  with  higher  smoke  and  fumes  from 
flame-retarded  materials  in  well  developed  fires  argues  against  this  mechanism  of  fire  resistance.  In 
the  following  sections  we  explore  the  molecular  basis  for  the  macroscopic  fire  parameters  in  an 
effort  to  make  quantitative  predictions  about  the  fire  resistance  of  a  polymer  from  its  chemical 
structure. 

3.  KINETICS. 

The  continuum  level  treatment  in  the  previous  section  shows  that  a  high  decomposition  temperature 
delays  ignition  (equation  3)  and  lowers  heat  release  rate  (equations  13  and  16)  by  increasing  the 
stored  heat  at  polymer  gasification.  Charring  is  a  solid  phase  mechanism  of  fire  resistance  which 
limits  the  amount  of  combustibles  which  can  be  generated  during  thermal  degradation  and  burning. 
In  this  section  we  provide  a  mechanism  for  fuel  generation  in  charring  and  nonchamng  polymers 
using  thermal  degradation  kinetics  and  show  that  the  pyrolysis  temperature  is  essentially  a  kinetic 
parameter. 

The  elementary  fuel  generation  step  is  thermal  degradation  of  the  polymer  [49,  51].  Typically,  it  is 
the  low  fraction  and  rate  of  production  of  volatile  fuel  at  fire  temperatures  and  the  low  heat  of 
combustion  of  this  fuel  which  make  polymers  intrinsically  fire  resistant.  Short-term  thermal 
stability  and  reduced  fuel  fraction  (increased  char  yield)  are  achieved  by  eliminating  hydrogen 
atoms  from  the  polymer  molecule  so  that  recombination  of  carbon  radicals  to  form  char  during 
thermal  degradation  is  kinetically  favored  over  hydrogen  abstraction/termination  reactions  which 
produce  volatile  fuel  fragments  [52,  53].  A  low  heat  of  combustion  of  the  volatile  thermal 
decomposition  products  is  achieved  by  substituting  heteroatoms  (e.g.,  halogens,  nitrogen, 
phosphorus,  sulfur,  silicon,  boron,  and  oxygen)  for  carbon  and  hydrogen  in  the  polymer  molecule. 
Heteroatoms  form  stable  gas  phase  combustion  products  which  are  either  low  in  fuel  value  (i.e.,  N2, 
SO2,  hydrogen  halides)  or  thermally  stable  solid  oxides  (i.e.,  Si02,  P2O5,  B2O3)  which  precipitate 
onto  the  polymer  surface  and  act  as  mass-  and  thermal-diffusion  barriers  [49, 50]. 

3.1  FUEL  GENERATION. 

A  description  of  the  rates  and  pathways  of  thermal  degradation  reactions  which  occur  in  the  solid 
state  during  burning  are  the  subject  of  this  section.  The  basic  thermal  degradation  mechanism 
leading  to  volatile  fuel  generation  in  polymers  has  been  described  [52, 54]  as  a  generalized  chemical 
bond  scission  process  consisting  of  primary  and  secondary  decomposition  events  as  illustrated 
schematically  in  figure  6.  The  primary  decomposition  step  can  be  main-,  end-,  or  side-chain 
scission  of  the  polymer.  Subsequent  thermal  degradation  reactions  depend  largely  on  the  chemical 
structure  of  the  polymer  [55, 56]  but  typically  proceed  by  hydrogen  fransfer  to  a-  or  (3-carbons, 
nitrogen  or  oxygen,  intramolecular  exchange  (cyclization),  side-chain  reactions,  small-molecule 
(SO,,  CO,,  S2)  elimination,  molecular  rearrangement,  and/or  unzipping  to  monomer.  Unzipping  or 
depolymerization  of  vinyl  polymers  is  characterized  by  a  kinetic  chain  length  or  “zip  length  which 
is  the  average  number  of  monomer  units  produced  by  a  decomposing  radical  before  the  radical  is 
deactivated  by  termination  or  transfer  [51].  Mathematically,  the  zip  length  is  the  ratio  of  the  rate 
constants  for  initiation  to  termination.  Aromatic  backbone  polymers  such  as  polycarbonate, 
polyimide,  and  polyphenyleneoxide  tend  to  decompose  in  varying  degrees  to  a  carbonaceous  char 
residue  through  a  complex  set  of  reactions  involving  crosslinking  and  bond  scission.  A  generally 
applicable,  detailed  mechanism  for  thermal  degradation  of  aromatic  polymers  is  unlikely. 

The  present  approach  avoids  the  need  for  a  detailed  degradation  mechanism  by  generalizing  thermal 
decomposition  as  a  transient  mass  balance  between  die  polymer,  a  reactive  intermediate,  and  the 


10 


degradation  products-fuel  gases  and  char.  Thus,  we  are  interested  in  the  products  of  thermal 
degradation  only  in  terms  of  their  contribution  to  the  heat  released  during  gas  phase  combustion. 
The  first  stage  of  thermal  degradation  produces,  in  a  generic  sense,  primary  volatiles  (gas  and  tar) 
and  possibly  a  primary  char  residue.  If  a  primary  char  forms,  further  decomposition  occurs  by 
dehydrogenation  to  form  the  secondary  gas  (principally  hydrogen)  and  a  thermally  stable 
secondary  carbonaceous  char. 


Polymer 


Reactive  intermediates 


lamt 


Volatiles 
(gas  +  tar) 


1 


f 


Primary 

ctjar 


Gas  Secondary 

(primarily  char 

hydrogen)  (primarily 

carbon) 


FIGURE  6.  GENERALIZED  THERMAL  DEGRADATION  MECHANISM  OF  POLYMERS 

Figure  7  clearly  shows  the  two-stage  process  of  thermal  decomposition  for  the  commercial 
polybenzimidazole  2,2,-(m-phenylene)-5,5/-bibenzimidazole  (CELAZOLE  PBI,  Hoechst  Celanese). 
Plotted  in  figure  7  is  the  mass  loss  rate  versus  temperature  curve  of  polybenzimidazole  at  a  linear 
heating  rate  of  10  K/min  in  an  inert  (nitrogen)  environment.  The  mass  loss  rate  peaks 
corresponding  to  the  primary  and  secondary  charring  reactions  are  evident. 


FIGURE  7.  THERMOGRAVIMETRIC  SCAN  OF  POLYBENZIMIDAZOLE  AT  A  LINEAR 
HEATING  RATE  OF  10  K/min  UNDER  NITROGEN 
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Decomposition  schemes  which  account  for  some  or  all  of  these  pyrolysis  products  (gas,  tar, 
primary  char,  secondary  char,  secondary  gas)  have  been  proposed  wherein  the  decomposition  steps 
occur  sequentially  (series),  simultaneously  (parallel),  or  in  some  combination  of  series/parallel  steps 
[37,  57,  58,  59,  60,  61,  62,  63,  64,  65,  66,  67].  Three  of  these  mechanistic  pyrolytic  reaction 
schemes  for  cellulose  have  been  reviewed  recently  [66]  including  the  single-step  first-order  model 
[65],  an  uncoupled  three-step  parallel  model  [60],  and  a  coupled  three-step  series-parallel  model 
[58, 67].  The  single-step  and  three-step  uncoupled  models  have  a  fixed  char  yield  as  an  adjustable 
parameter  while  the  three-step  coupled  model  has  a  variable  char  yield.  All  of  the  models  predict 
rate-dependent  peak  decomposition  temperatures.  Variable  (n-th)-  order  decomposition  kinetics 
have  been  fit  to  mass  loss  data  for  char-forming  polymers  with  reasonable  success  using  reaction 
order  [68]  and  empirical  weighting  factors  [69]  as  adjustable  parameters.  However,  little  insight  is 
gained  into  the  reaction  pathways  from  these  curve-fitting  exercises. 

A  simple,  solid-state  fuel  generation  model  has  been  derived  [70]  from  the  following  assumptions 
about  the  process  of  polymer  thermal  degradation  in  fires: 

1.  A  reactive  intermediate  I*  is  generated  in  the  polymer  dissociation  (initiation)  step  which  is 
in  rapid  dynamic  equilibrium  with  the  parent  polymer,  P,  but  is  consumed  in  the  process  of 
gas  and  char  formation  such  that  its  concentration  never  becomes  appreciable  and  decreases 
slowly  over  time  as  the  polymer  is  consumed.  Consequently,  the  rate  of  change  of  reactive 
intermediate  with  time  is  insignificant  compared  to  the  rate  of  polymer  consumption,  gas 
production,  and  char  formation  so  that  for  computational  purposes  the  change  in 
intermediate  concentration  with  time  is  neglected.  This  is  the  stationary-state  hypothesis. 

2.  Thermal  decomposition  of  primary  char  to  secondary  char  and  gas  is  slow  compared  to  the 
formation  of  the  primary  char  at  typical  flaming  surface  temperatures  of  400-700°C. 
Consequently,  only  the  primary  char  is  considered  in  formulating  the  reaction  set  for  the 
mass  loss  model. 

3.  The  thermal  degradation  environment  in  the  pyrolysis  zone  of  a  burning  solid  polymer  is 
non-oxidizing  or  anaerobic.  Dissolved  molecular  oxygen  and  oxygen  diffusion  into  the 
pyrolysis  zone  of  the  solid  are  considered  negligible  with  respect  to  their  effects  on  gas  and 
char  formation  so  that  solid-state  oxidation  reactions  can  be  neglected  in  the  fuel  generation 
model  for  polymers  in  fires.  This  assumption  does  not  preclude  the  possibility  of  surface 
mass  loss  due  to  thermoxidative  reactions  at  the  polymer-air  interface  under  nonflaming 
conditions,  e.g.,  thermogravimetric  (TGA)  experiments  conducted  in  air  or  smoldenng 
combustion. 

Figure  8  shows  data  [71]  for  a  variety  of  pure,  unfilled  polymers  plotted  as  the  char  yield  measured 
after  flaming  combustion  in  a  fire  calorimeter  versus  the  char  residue  at  900  ±  100°C  for  the  same 
material  after  anaerobic  pyrolysis.  It  is  seen  that  the  char  yield  of  a  material  in  a  fire  is  essentially 
equal  to  its  residual  mass  fraction  after  pyrolysis  in  an  oxygen-free  environment  at  temperatures 
representative  of  the  char  temperature  in  a  fire.  Although  oxidative  degradation  products  have  been 
identified  at  the  surface  of  noncharring  olefinic  polymers  after  flaming  combustion  [56,  72],  the 
data  in  figure  8  suggest  that  oxidation  reactions  are  insignificant  in  the  pyrolysis  zone  of  a  burning 
polymer  as  evidenced  by  the  close  agreement  between  fire  char  yield  and  anaerobic  pyrolysis 
residue. 
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FIGURE  8.  FIRE  CHAR  YIELD  VERSUS  ANAEROBIC  PYROLYSIS  RESIDUE  FOR  A 

VARIETY  OF  POLYMERS 


The  generalized  combustion  and  pyrolysis  schemes  of  figures  1  and  6,  respectively,  in  combination 
with  assumptions  1-3  lead  to  the  simplified  kinetic  model  for  polymer  burning  which  is  shown  in 
figure  9. 
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FIGURE  9.  THERMAL  DEGRADATION  MODEL  FOR  POLYMERS  IN  FIRES 

This  simplified  scheme  reduces  thermal  degradation  of  polymer  P  to  a  single  step  involving  parallel 
reactions  of  an  active  intermediate  I*  to  form  to  gas  G  and  char  C.  In  figure  9,  kj  is  the  rate 
constant  for  initiation,  and  kr,  kg,  and  kc  are  the  rate  constants  for  termination  by  recombination  (kr), 
hydrogen  transfer  to  gaseous  species  (kg),  and  crosslinking  to  char  (kc),  respectively.  The  rate 
constants  are  assumed  to  have  an  Arrhenius  temperature  dependence  of  die  form  k(T)  =  A 
exp[-E/RT]  where  A  and  E  are  the  frequency  factor  and  activation  energy,  respectively,  at 
temperature  T.  Neglecting  solid-state  oxidation  the  thermal  decomposition  reactions  are 
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( rapid  equilibrium )  (18) 

(slow)  (19) 

(slow)  (20) 

and  the  system  of  rate  equations  for  the  species  at  time,  t,  is 


f = -k,p  +V 

(21) 

^  =  k,P-(k,  +  k,+  lc)l» 

(22) 

dG  _  t  t* 

dt  ~  s 

(23) 

dC  _  T^i 
dt  '  ^ 

(24) 

According  to  the  stationary-state  hypothesis,  dl*/dt  «  0  so  that  equation  22  provides  the  useful 
result 


I* 


ki 

kr  +  kg  +  kc 


P 


KP 


where  K  =  k;  AX  +  kg  +  kc)  is  the  pseudo-equilibrium  constant  for  the  polymer  dissociation 
reaction.  As  the  ratio  of  initiation  to  termination  rate  constants,  K  represents  the  kinetic  chain  length 
for  degradation  by  depolymerization  [51, 73].  Substituting  I*  =  K  P  into  equations  21, 23,  and  24, 

^  =  -[k,-Kk,]P  (25) 

^  =  k.KP  (26) 

at  6 

(27) 

(28) 


f  -k*p 


With  I*  «  P,  G,  C,  the  total  mass  balance  in  terms  of  the  initial  mass,  mo,  is 


mr 


=  P  +G  +  C+I*  =  P+G  +  C 
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From  equations  25-28  with  dm0/dt  =  0 


dP 

dt 


dC_dG 
dt  dt 


=  -[k,-Kkr]P 


(29) 


The  sensible  mass  of  the  sample  as  measured,  for  example,  in  a  TGA  experiment  or  fire  calorimeter 
testis 


and  with  equation  26 


m  =  P  +  C  + 1*  *  P  +  C 
dm  dP  dC  dG 


dt 


dt  +  dt 


dt  -  -Kk,P 


(30) 


Equation  29  can  be  solved  immediately  for  P  in  the  isothermal  case  with  initial  condition  P  =  P0  = 

mo  @  t  =  0, 


P  =  m0exp|-[ki-Kkr]t| 

Substituting  the  isothermal  result  for  P  into  equation  30  and  separating  variables 

J  dm'  =  -J^  Kkgm0exp(-kpt)dt 


(31) 


(32) 


where  kp  in  the  exponential  of  the  integrand  on  the  right-hand  side  of  equation  32  is  the  overall  rate 
constant  for  pyrolysis  and  is  assumed  to  have  the  Arrhenius  form 


kp  =  kj-Kk,.  =  K(kg  +  kc)  =  Aexp 


V 

RT 


(33) 


in  terms  of  the  global  activation  energy  Ea  and  frequency  factor  A  for  pyrolysis.  The  isothermal 
solution  of  equation  32  is 


or 


m(t) 

mn 


m(t) 

m„ 


=  1- 


kg  +  kc 


(l  -  e-kpt ) 


=  Yc(D  +  [l-  YefDje^p* 


(34a) 


(34b) 


Equation  34  shows  that  as  t  -»  °°  the  residual  mass  approaches  an  equilibrium  value  at  constant 
temperature  given  by 


YC(T)  = 


m(°°) 

m„ 


kg  +  kc 


(35) 
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where  YC(T)  is  the  equilibrium  residual  mass  fraction  or  char  yield  at  temperature  T  in  terms  of  the 
rate  constants  for  gas  and  char  formation.  Equation  35  predicts  a  finite  char  yield  at  infinite  time  if 
kc  >  0  and  zero  char  if  kc  =  0. 

The  physical  significance  of  a  temperature-dependent,  equilibrium  char  yield  as  the  ratio  of  rate 
constants  for  gas  and  char  formation  is  consistent  with  the  use  of  group  contributions  for  the  char 
forming  tendency  of  polymers  developed  by  Van  Krevelen  [52,  53]  (see  the  following  section).  If 
kg  and  kc  have  Arrhenius  forms,  equation  35  can  be  written 


YC(T)  = 


1  +  — 1  exp  f-(Eg-Ec)/RT 
A.  o 


(36) 


where  Ec,  Eg  and  Ac,  Ag  are  the  activation  energies  and  frequency  factors  for  char  and  gas 
formation,  respectively.  The  crossover  temperature,  Tcr,  is  defined  [66]  as  the  temperature  at  which 
the  rates  of  gasification  and  crosslinking  are  equal,  i.e.,  when  kg  =  k^ 


.  (Eg-Ec) 

cr  Rln[Ag/Ac] 


(37) 


It  follows  from  equation  35  that  the  crossover  condition,  kg  =  kc,  corresponds  to  the  equilibrium 
residual  mass  fraction,  Yc(Tcr)  =  0,50. 

If  YC(T)  is  the  char  yield  at  a  temperature  above  the  major  mass  loss  transition  or  is  independent  of 
temperature  then  YC(T)  =  |i  =  constant  and  equation  34  is  the  solution  for  the  isothermal  mass  loss 
history  of  a  filled  polymer  with  a  nonvolatile  mass  fraction  p  satisfying  the  rate  law 


^  =  -kp(m-|im0)  (38) 

although  equation  38  was  not  assumed  a  priori  in  the  present  derivation. 

The  previous  results  were  derived  for  the  isothermal  (constant  temperature)  case  but  many 
processes  of  interest  are  conducted  under  nonisothermal  conditions,  e.g.,  thermogravimetric 
analyses  at  constant  heating  rate  or  fuel  generation  in  the  pyrolysis  zone  of  a  burning  polymer.  To 
calculate  the  instantaneous  mass  fraction  m(t)/m0  during  a  constant  heating  rate  experiment  where 
dT/dt  =  constant  =  P,  begin  by  eliminating  P  between  equations  29  and  31  and  integrating 

f  dm'  =  (1-YJ  f  dP'  (39) 

J  m0  J  Pq 


or  since  P0  =  m0, 


m(T) 

m0 


YC(T)  +  [1  -  YC(T)] 


PCD 

Po 


For  nonisothermal  conditions  P(T)/P0  in  equation  40  is  obtained  from  equation  29  as 


(40) 
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where  the  constant  heating  rate  P  =  dT/dt  transforms  the  variable  of  integration  from  time  t  to 
temperature  T;  and  A,  Ea  are  the  global  frequency  factor  and  activation  energy  of  pyrolysis, 
respectively. 

The  right-hand  side  of  equation  41  is  the  exponential  integral  which  has  no  closed  form  solution. 
An  approximate  solution  for  the  exponential  integral  which  is  accurate  to  within  2%  over  the  range 
of  activation  energies  and  temperatures  encountered  in  thermal  analysis  and  combustion  is  [74] 


-ART2  [  E,| 
P  (Ea  +  2RT)  6Xp  \  RTj 


-kpRT2 
P(Ea  +  2RT) 


(42) 


Defining 

kpRT2 

y  P(Ea  +  2RT) 


(43) 


the  solution  of  equation  41  takes  the  form 


(44) 


Substituting  equation  44  into  equation  40,  the  residual  mass  fraction  in  a  constant  heating  rate 
experiment  is 


m(T) 

m0 


=  YC(T)  +  [1  -  Yc(T)]e-y 


(45) 


which  is  the  same  form  as  the  isothermal  solution  (equation  34).  Recent  work  [75]  on  the  thermal 
degradation  of  polymers  in  fires  shows  that  pseudo  first-order  mass  loss  kinetics  (e.g.,  equation  45) 
are  a  good  approximation  to  n-th  order  kinetics  at  constant  heating  rate  for  typical  reaction  orders 
n  <  4. 


Figure  10  shows  the  agreement  of  equation  45  with  nonisothermal  experimental  data  for  phenolic 
triazine,  a  char-forming  thermoset  polymer  with  Ea  =  178  kJ/mol  and  A  =  109  sec-1  determined 
from  isothermal  weight  loss  experiments.  Solid  lines  in  figure  10  are  the  fits  of  equation  45  to 
experimental  data  (open  symbols)  at  linear  heating  rates  P  =  1, 5,  and  20  K/min.  The  dashed  line  is 
die  temperature  dependent  residual  mass  YC(T)  calculated  from  equation  36  for  an  energy  barrier  to 
char  formation  E£  -  Ec =  +30  kJ/mol  and  frequency  factor  ratio  AJAC=  17  determined  in  separate 
isothermal  experiments  [70].  The  agreement  of  this  simple  model  with  experimental  data  is  seen  to 
be  very  good  for  this  polymer  up  to  temperatures  =700°C  where  secondary  char  formation,  which 
was  neglected  in  the  model,  begins. 
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FIGURE  10.  COMPARISON  OF  THERMAL  DEGRADATION  MODEL  TO 
NONISOTHERMAL  TGA  DATA  (NITROGEN  PURGE)  FOR  A  PHENOLIC  TRIAZINE 

THERMOSET  POLYMER 


The  fractional  mass  loss  rate  during  a  linear  temperature  ramp  is  obtained  by  differentiating 
equation  45  with  respect  to  time, 


- 1  <im(T) 
m0  dt 


+ 


(46) 


=  (l  -  YC(T))  kp(T)  e  -y  +  p  YC(T)  (l  -  YC(T))  ^^(l  - e  "y) 


Since  the  rate  of  change  of  YC(T)  is  small  compared  to  the  fractional  mass  loss  rate  at  pyrolysis  (see 
figure  10),  the  approximation  YC(T)  =  p  =  constant  is  made  so  that  dYc  /dt  =  0  and  equation  46 

simplifies  to 


- 1  dm(T) 
m0  dt 


(l-p)kpe  y 


(47) 


The  peak  value  of  the  fractional  mass  loss  rate  can  be  found  by  differentiating  equation  47  with 
respect  to  time  and  setting  this  second  derivative  of  the  residual  mass  fraction  equal  to  zero, 


d2  j 

m(T)\ 

dt2 ' 

mo  j 

=  P(l-H) 


d  T 


kpe 


i  =  M 


kpe 


-y 


PEa 


RT 


=  0 


(48) 


Equation  48  has  two  roots:  the  trivial  case  p  =  1  and 


kp(max) 


(49) 
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where  Tp  is  the  temperature  at  maximum  mass  loss  rate  during  the  course  of  the  linear  heating 
history.  For  example,  Tp  =  635°C  at  the  peak  fractional  mass  loss  rate  of  0.23  mg/g-s  for  the 
primary  decomposition  of  polybenzimidazole  at  a  constant  heating  rate  of  10  K/min  in  figure  7. 


An  analytic  result  for  the  peak  fractional  mass  loss  rate  in  a  constant  heating  rate  experiment  is 
obtained  by  substituting  equation  49  into  equation  47 


- 1  dm 


m0  dt 


max 


=  t1'!*) 


PE. 


e'RT' 


(50) 


where  the  exponent  r  of  the  natural  number  e  in  the  denominator  has  the  value 


r  = 


1  +  2R^ 


1-1 


(51) 


For  the  usual  case  where  Ea »  2RTp,  equation  50  simplifies  to 


- 1  dm 
m0  dt 


max 


-  (M 


PEa 

eRT„2 


(52) 


Once  A  and  Ea  have  been  determined,  the  heating  rate  dependent  temperature  at  peak  mass  loss  rate 
Tp  is  obtained  from  the  root  Ea/RTp  of  equation  49  written  in  the  form 


In 


Ea 

2 

l 

Ea' 

|  1 « 

PR 

RTp 

T 

.3 

t  in 

AEa 

0 


(53) 


The  effect  of  heating  rate  on  peak  mass  loss  temperature  Tp  embodied  in  equation  53  is  not 
insignificant  when  comparing  laboratory  ((3  ~  10  K/min)  and  fire  ((3  ~  300  K/min)  heating  rates. 
To  illustrate  this  effect,  pyrolysis  experiments  were  conducted  under  a  nitrogen  atmosphere  in  a 
commercial  thermogravimetric  analyzer  at  various  linear  heating  rates  [70].  Figure  11  shows  the 
experimental  data  for  the  temperature  at  peak  mass  loss  rate  Tp  for  samples  weighing 
<  5  milligrams.  Polymers  tested  include  a  linear  semicrystalline  thermoplastic  polyethylene  (PE), 
Ea  =  264  kJ/mol,  p  =  0  [76,  77],  which  degrades  to  short  hydrocarbon  chains  by  random  chain 
scission;  a  Unear  amorphous  thermoplastic  polymer  poly(methylmethacrylate)  (PMMA), 
Ea  =  160  kJ/mol,  p  =  0  [76,  77],  which  thermally  degrades  to  monomer  by  an  unzipping 
mechanism;  and  a  highly  crosslinked,  amorphous,  thermoset  polymer  phenoUc  triazine  (PT),  Ea  = 
178  kJ/mol,  p  =  0.7  [70, 74],  which  degrades  by  decyclization,  crosslinking,  and  charring. 

The  accuracy  of  equation  52  was  tested  by  measuring  the  peak  mass  loss  rate  at  the  decomposition 
temperatures  plotted  in  figure  11  for  PMMA,  PE,  and  PT.  The  measured  peak  mass  loss  rates  are 
compared  in  figure  12  to  the  peak  mass  loss  rates  calculated  using  equation  52  with  the  previous 
activation  energies  and  pyrolysis  residues  for  these  polymers  along  with  the  peak  mass  loss 
temperatures  plotted  in  figure  11.  Close  proximity  of  the  data  to  the  equivalence  line  on  the  log-log 
plot  indicates  good  agreement  between  measured  and  calculated  peak  mass  loss  rate  using 
equation  52  at  heating  rates  ranging  from  (3  =  1  to  200  K/min. 
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FIGURE  11.  PEAK  MASS  LOSS  TEMPERATURE  VERSUS  HEATING  RATE  FOR 
POLYMETHYLMETHACRYLATE  (PMMA),  POLYETHYLENE  (PE),  AND 

PHENOLIC  TRIAZINE  (PT) 


FIGURE  12.  CALCULATED  VERSUS  MEASURED  PEAK  MASS  LOSS  RATE  IN  TGA 
FOR  PMMA,  POLYETHYLENE,  AND  PHENOLIC  TRIAZINE  AT  INDICATED  HEATING 
RATES,  |3  =  1, 3, 5, 10, 20, 30, 100,  and  200  K/min 

3.2  DECOMPOSITION  TEMPERATURE. 

Intrinsic  fire  resistance  requires  short-term  thermal  stability  at  fire  temperatures  as  discussed  in 
previous  sections.  Van  Krevelen  [52, 73]  identified  a  number  of  experimental  indices  which  can  be 
used  to  characterize  the  short-term  thermal  stability  of  a  polymer.  The  indices  which  are  determined 
from  temperature  scanning  thermogravimetry  experiments  at  linear  heating  rates  P  ~  3-10  K/min 
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are  Td;0,  the  temperature  in  Kelvin  at  which  the  polymer  weight  loss  is  just  measurable;  Td,i/2,  the 
temperature  at  which  the  weight  loss  reaches  50%  of  its  final  value;  Td>max,  the  temperature  at  the 
maximum  rate  of  weight  loss;  and  E*  the  activation  energy  for  pyrolysis.  Noting  that  all  of  these 
temperatures  depend  on  heating  rate  and  that  Td)inax  is  identical  to  Tp  in  the  present  work,  Van 
Krevelen  proposed  the  following  interrelationships  based  on  tabulated  data  for  37  polymers  (all 
temperatures  are  Kelvin): 


T  =  0.9  T 

d,0  d.1/2 


(54) 


T  =  T  ~  T 

d,  max  p  d,l/2 


(55) 


Tdm  *  423  +  Ea  (kJ/mol) 


(56) 


Equations  54-56  imply  that  T  or  Td  m  can  be  used  as  the  primary  indicator  of  thermal  stability  with 
the  other  indices  being  roughly  estimated  from  it.  A  plot  of  Tp  versus  Ea  from  thermogravimetric 
analyses  of  various  polymers  [52,  73]  at  heating  rates  on  the  order  of  (3  =  3  K/min  is  shown  in 
figure  13.  The  trend  line  is  in  qualitative  agreement  with  equation  56,  but  the  correlation  coefficient 
is  too  low  (<  0.5)  for  predictive  purposes. 


Activation  Energy  for  Pyrolysis,  Ea  (kJ/mol) 


FIGURE  13.  DEPENDENCE  OF  DECOMPOSITION  TEMPERATURE  ON  ACTIVATION 

ENERGY  FOR  PYROLYSIS 


Despite  the  general  trend  of  Tp  increasing  with  Ea  in  figure  13,  the  scatter  is  relatively  high  because 
of  the  variation  in  the  frequency  factor  A  between  materials  and  its  influence  on  Ea/Tp  embodied  in 
equation  53.  In  the  following  paragraphs  we  derive  the  proportionality  between  Tp  and  Ea  from 
Arrhenius  kinetics  (i.e.,  equation  53)  and  show  that  it  is  material  specific. 

The  Arrhenius  relation  for  the  temperature  dependence  of  reaction  rate  constants,  e.g.,  kp  =  A 
exp[-Ea/RT],  was  originally  derived  with  the  assumption  that  the  reactants  are  at  equilibrium  with  a 
high-energy  intermediate  state  and  could  proceed  to  products  with  no  further  energy  requirements 
[78,  79].  Dynamic  equilibrium  between  the  dissociating  polymer  P  and  activated  complex  I* 
(equation  18)  satisfies  this  condition  and  implies  a  thermodynamic  significance  for  the  pyrolysis 
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kinetic  parameters  A  and  Ea.  From  the  definition  of  the  free  energy  at  equilibrium  and  equation  33, 


In  K  =  In 


kp 

1 

1 

* 

O 

<1 

l 

AH* 

AS* 

kg  +  kc 

RT 

RT 

R 

(57) 


where  as  previously,  kp  is  the  global  rate  constant  for  pyrolysis  and  kg  and  kc  are  the  rate  constants 

for  gas  and  char  formation,  respectively.  In  equation  57,  AG*,  AH*,  and  AS  are  the  molar  free 
energy,  enthalpy,  and  entropy  of  pyrolysis,  respectively.  Assuming  Arrhenius  forms  for  kp,  kg,  and 
k,.,  it  follows  that 

kp  =  Ae"^7  = 


(kg  +  ^c)  * 


AS  /R 


-AH  /RT 


(58) 


where  A  =  (kg  +  kc)e  AS*/R  =  (kp/K)  e  AS*/R  is  the  pre-exponential  factor,  and  Ea  =  AH*  is  the 
global  activation  enthalpy  for  anaerobic  mass  loss  (pyrolysis).  From  equation  58  with  kp  = 
(3Ea/RTp2  (equation  49)  at  peak  mass  loss  rate,  the  entropy  change  for  the  solid  ->  gas  phase 
transition  at  temperature  Tp  is 


AS*(P)  =  Rln 


(kart'/e^ 

p  I 


(59) 


If  AG*  =  0  for  this  [dynamic]  equilibrium  process  then  K  =  1  and  a  molecular  heating  rate  can  be 
defined  having  units  (K/s) 


The  molecular  heating  rate  (3*  is  a  constant  if  changes  in  A,  Ea,  or  their  ratio  A/Ea  compensate  [80] 
for  changes  in  Tp2  with  p  (see  figure  11).  Rewriting  equation  53  as 


Ea 

RTp 


In 


artp! 

PE. 


and  substituting  P* 
be  expressed  as 


=  ARTp2/Ea  into  the  logarithmic  term,  the  peak  pyrolysis  rate  temperature  can 


E» 

R  ln[p*/p] 


(60a) 


Figure  14  shows  the  data  of  figure  11  plotted  as  Tp  versus  l/ln(P*/p)  according  to  equation  60a 
with  p*  adjusted  to  obtain  slope  Ea/R  for  PE  (Ea  =  264  kJ/mol),  PT  (Ea  =  178  kJ/mol),  and  PMMA 
(E  =  160  kJ/mol).  It  is  seen  that,  to  within  the  expected  accuracy  of  Ea/R  (±  10%),  the  intercept  in 
figure  14  which  is  unconstrained  in  the  fitting  exercise  to  obtain  P*  passes  through  the  origin  as 
required  by  equation  60a.  The  correlation  coefficient  for  all  of  the  plots  in  figure  14  is  better  than 
99%  over  the  relatively  narrow  range  of  heating  rates  P  =  1-200  K/min  fitted  for  the  extrapolation. 
From  the  P*  in  figure  14  and  typical  RT  2/Ea  =  20  K,  the  frequency  factors  for  pyrolysis  are 
A  =  P*/20  K  =  1017  -  1012  s"1,  typical  of  frequency  factors  encountered  in  studies  of  polymer 
pyrolysis  [81,82]. 
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FIGURE  14.  PEAK  PYROLYSIS  TEMPERATURE  VERSUS  RECIPROCAL  OF  THE 
NATURAL  LOGARITHM  OF  THE  REDUCED  HEATING  RATE  ACCORDING  TO 
EQUATION  60  (Dashed  lines  have  slope  Ea/R  for  indicated  (5*.) 

Substituting  [3*  =  ARTp2/Ea  into  equation  59  gives  AS(P)  =  R  ln[(3*/(3],  which  has  the  same  form 
as  the  statistical  thermodynamic  entropy  if  P*/p  is  the  equilibrium  number  of  states  having  energy 
Ea  at  temperature  Tp.  The  analogy  between  (3*/(3  and  N  can  be  shown  by  defining  temperatures  in 
the  vicinity  of  the  peak  pyrolysis  temperature  T,  =  (Tp  -  AT/2)  and  T2  =  (Tp  +  AT/2)  with  AT  the 
width  of  die  pyrolysis  peak  such  that  AT  «  Tp.  It  is  found  (c.f.,  method  of  section  5)  that  for  a 
first-order  decomposition  process  at  constant  heating  rate  the  characteristic  width  of  the  mass  loss 
transition  (pyrolysis  peak)  in  degrees  Kelvin  is  AT  ~  RTp2/Ea.  Thus  P*  =  ARTp2/Ea  =  AAT.  The 
time  interval  corresponding  to  this  temperature  range  is  At  =  AT/p  so  that  P*/p  «  (AAT)/(AT/At)  = 
AAt  =  N  where  N  is  the  number  of  molecular  collisions  occurring  in  time  At  at  temperatures  in  the 
vicinity  of  Tp  with  sufficient  energy  Ea  to  spontaneously  produce  gas  and/or  char. 

Thus,  AS(P)  =  R  ln[P*/p]  =  R  ln[N]  is  consistent  with  the  statistical  thermodynamic  interpretation 
of  entropy  as  the  logarithm  of  the  number  of  excited  states  at  Tp  and  equation  60a  can  be  written  in 
the  various  forms 


Ea  _  AH*  AH* 

Tp(P)  "  R  ln[p*/p]  “  Rln[N]  “  AS*(P) 


(60b) 


The  inverse  proportionality  between  decomposition  temperature  Tp  and  ln[N]  is  physically  realistic 
because  the  greater  the  number  of  collisions  resulting  in  mass  loss  (low  heating  rate  =  large  At  = 
large  ln[N])  the  earlier  the  peak  mass  loss  temperature  will  occur  in  the  heating  program 
(low  Tp)  since  the  fractional  mass  loss  for  complete  pyrolysis  is  independent  of  heating  rate 
(see  section  4.1). 

Although  equation  60a  is  easily  recast  in  the  form  of  equation  56  by  defining  a  reference 
temperature  T0  =  423  K  corresponding  to  a  reference  heating  rate  P  =  P0,  it  is  not  sufficiently 
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accurate  for  predictive  purposes.  For  this  reason.  Van  Krevelen  extended  the  use  of  additive  molar 
group  contributions  which  he  had  previously  used  to  calculate  the  equilibrium  glass  transition  and 
melting  temperatures  to  the  calculation  of  the  thermal  decomposition  temperature  [52,  91].  Additive 
molar  group  contributions  had  long  been  used  for  quantitative  prediction  of  the  equilibrium 
properties  of  ordinary  molecules  and  polymers,  e.g.,  heats  of  combustion  [83],  thermophysical  and 
transport  properties  [84,  85],  chemical  reactivity  [86],  and  biological  activity  [87].  In  Van 
Krevelen’ s  method  of  calculating  the  thermal  decomposition  temperature,  the  molar  thermal 
decomposition  function  Yd  m  is  related  to  the  product  of  the  polymer  repeat  unit  molecular  weight 
M  and  the  decomposition  temperature  viz. 

M-Td  1/2  =  Yd  1/2  (61) 

From  the  previous  discussion  it  is  apparent  that  Td  1/2  =  Tp  is  a  kinetic  (rate-dependent)  parameter 
and  so  the  molar  thermal  decomposition  function  must  be  also.  Comparison  of  equations  60b  and 
61  show  that  for  a  particular  heating  rate  (30  the  molar  decomposition  function  is  related  to  the  molar 
heat  AH*  and  entropy  AS*((30)  of  decomposition  as 


=  Yd,1/2(P0) 


M-AH* 

AS*(p0) 


The  decomposition  temperature  is  calculated  by  summing  the  individual  molar  decomposition 
functions  Yd  m  for  each  chemical  group  of  the  molar  mass  Mj  comprising  the  repeat  unit 


Td,  1/2  “ 


Yd,  1/2 

M 


i^l  ni  Yd,  1/2 


.f^njMi 


(62) 


Equation  62  is  the  same  form  as  those  used  to  predict  the  glass  transition  and  melting  temperatures 
of  polymers  using  the  appropriate  values  for  the  molar  group  contributions  [52].  The  empirical 
method  of  Van  Krevelen  for  calculating  thermal  decomposition  temperature  requires  a  molar  group 
contribution  yJ  m  for  each  of  the  chemical  groups  comprising  the  polymer.  These  molar  group 
contributions  were  determined  for  decomposition  temperatures  measured  at  a  linear  heating  rate 
P  «  3  K/min.  The  following  example  is  illustrative  of  the  additive  method  of  group  contributions 

developed  by  Van  Krevelen. 

To  calculate  the  decomposition  temperature  of  the  fire-resistant  commercial  polymer  CELAZOLE 
PBI  (2,2'-(m-phenylene)-5,5'-bibenzimidazole),  first  write  the  chemical  repeat  unit  structure  for  PBI 


Comparing  the  chemical  structure  to  tabulated  group  contributions  for  Yd  j/2  >  the  data  in  table  1 
can  be  assembled.  From  table  1  and  equation  62  the  peak  decomposition  temperature  of  PBI  at  low 
heating  rate  is 
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(2)(105)  +  (1-65) 
(2)(116.12) +  76.09 


•  1000 


892  K 


which  is  in  general  agreement  with  the  literature  value  [52]  Tp  «  903  K  and  the  primary 
decomposition  peak  in  figure  7. 

TABLE  1.  GROUP  CONTRIBUTIONS  TO  THE  DECOMPOSITION  TEMPERATURE  OF 
POLYBENZIMIDAZOLE  AFTER  VAN  KREVELEN  [91] 


Group 

Ouantitv 

M] 

(g/mol) 

Yd,  1/2 
(K-kg/mol) 

2 

116.12 

105 

1 

76.09 

65 

While  calculation  of  Tp  for  commercial  PBI  is  straightforward,  it  is  often  the  case  that  the  molar 
group  contributions  are  not  available  in  the  literature  for  predicting  decomposition  temperatures  of 
new,  thermally  stable  polymers  containing  novel  backbone  and  pendant  structures.  Recognizing 
this  limitation  of  the  group  contribution  method,  Bicerano  [84]  generalized  the  additive  scheme 
using  graph  theory  to  develop  atomic  connectivity  indices  which  replace  the  larger  chemical  groups 
in  the  traditional  additive  approach  as  the  principle  descriptors  of  the  topology  of  the  polymer  repeat 
unit.  Bicerano’ s  approach  is  still  empirical  in  that  the  connectivity  indices  must  be  determined  from 
experimental  data  or  by  correlation  with  known  group  contributions.  However,  once  the 
connectivity  indices  are  determined  for  a  particular  property,  that  property  may  be  predicted  for  any 
chemical  structure  for  which  the  atomic  (as  opposed  to  group)  indices  are  known.  The  connectivity 
index  method  of  predicting  polymer  properties  from  atomic  composition  can  be  accomplished  with 
a  hand  calculator  but  is  more  typically  implemented  using  a  stand-alone  computer  code 
(SYNTHIA)  or  accessed  as  a  module  in  molecular  modeling  software  (Molecular  Simulations,  Inc., 
San  Diego,  CA). 

Molecular  dynamics  simulations  of  polymer  thermal  degradation  have  also  recently  become 
available  [88].  Molecular  dynamics  simulations  are  the  most  fundamental  and  powerful  approach 
to  understanding  and  predicting  polymer  thermal  stability  because,  unlike  the  additive  approaches, 
molecular  dynamics  can  predict  the  thermal  decomposition  products  of  polymers  [89]. 

Table  2  is  a  listing  of  the  temperature  at  peak  mass  loss  rate  T  [90, 91]  and  the  pyrolysis  residue 
[X  [90, 91, 92]  at  850  ±  50°C,  for  the  indicated  polymers  heated  in  an  inert  environment  at  a  linear 
heating  rate  (3  =  3-10  K/min  according  to  standard  methods  [93].  Also  listed  are  the  limiting 
oxygen  index  [56,  90,  94]  and  the  Underwriters  Laboratory  UL  94  Vertical  or  Horizontal  bum 
ranking  [95].  The  limiting  oxygen  index  [96]  is  the  minimum  concentration  (volume  percent)  of 
oxygen  in  the  environment  required  to  sustain  candlelike  burning  (downward  flame  spread)  after 
ignition  of  a  3-  x  6-  x  150-mm  specimen  at  room  temperature  using  a  small  methane  flame.  A 
limiting  oxygen  index  above  20  would  indicate  self-extinguishing  behavior  when  ignited  at  ambient 
(20°C,  20  volume  percent  oxygen)  conditions. 


25 


TABLE  2.  PEAK  PYROLYSIS  TEMPERATURE  (Tp),  PYROLYSIS  RESIDUE  ([i),  LIMITING 
OXYGEN  INDEX  (LOI),  AND  UL  94  RANKING  OF  SOME  POLYMERS  (Values  in 

parentheses  are  calculated  estimates.) 


Polymer 

T, 

(°6 

(%) 

LOI 

(%) 

UL  94 
Ranking 

Polybenzobisoxazole  (PBO) 

789 

75 

56 

V-0 

Polyparaphenylene 

652 

75 

55 

(V-0) 

Polybenzimidazole  (PBI) 

630 

70 

42 

V-0 

Polyamideimide  (PAI) _ 

628 

55 

45 

V-0 

Polyaramide  (KEVLAR) 

628 

43 

28 

V-0 

Polyetherketoneketone  (PEKK) 

619 

62 

40 

V-0 

Polyetherketone  (PEK) 

614 

56 

40 

V-0 

Polvtetrafluoroethylene  (PTFE) 

612 

0 

95 

V-0 

Polyetheretherketone  (PEEK) 

606 

50 

35 

V-0 

Polyphenylsulfone  (PPSF) 

606 

44 

38 

V-0 

Polypara(benzoyl)phenylene  (PX) 

602 

66 

41 

V-0 

Fluorinated  Cyanate  Ester 

583 

44 

(40) 

V-0 

Polyphenylenesulfide  (PPS) 

578 

45 

44 

V-0 

Polyetherimide  (PEI) 

575 

52 

47 

V-0 

Polypyromellitimide  (PI) 

567 

70 

37 

V-0 

Liquid  Crystal  Polyester 

564 

38 

40 

V-0 

Polycarbonate  (PC) 

546 

25 

26 

V-2 

Polysulfone  (PSF) 

537 

30 

30  ' 

V-l 

Polyethylene  (PE) 

505 

0 

18 

HB 

Polyamide  6  (PA6) 

497 

1 

21 

HB 

Polyethylenenapthalate  (PEN) 

495 

24 

32 

V-2 

Polyphthalamide 

488 

3 

(22) 

HB 

Phenolic  Triazine  Cyanate  Ester  (PT) 

480 

62 

30 

V-0 

Polyethyleneterephthalate  (PET) 

474 

13 

21 

HB 

Cyanate  ester  of  Bisphenol-A  (BCE) 

470 

33 

24 

V-l 

Polydimethylsiloxane  (PDMS) 

444 

0 

30 

HB 

Acrylonitrile-butadiene-styrene  (ABS) 

444 

0 

18 

HB 

Polyurethane  elastomer  (PU) 

422 

3 

17 

HB 

Polymethylmethacrylate  (PMMA) 

398 

2 

17 

HB 

Polychlorotrifluoroethylene 

380 

0 

95 

V-0 

Polystyrene  (PS) 

364 

0 

18 

HB 

Polyoxymethylene  (POM) 

361 

0 

15 

HB 

Poly(a-methylstyrene) 

341 

0 

18 

HB 

Polyvinylchloride  (PVC) 

270 

11 

50 

V-0 

The  UL  94  test  [97]  ranks  materials  with  respect  to  their  time  to  self-extinguish  after  a  10-second 
Bunsen  burner  ignition.  A  V-0  rating  in  this  table  means  that  a  thin  (<  1-mm-thick)  specimen  of  the 
polymer  self-extinguished  in  a  vertical  orientation  within  10  seconds  after  removal  of  the  flame 
while  V-l  and  V-2  rankings  indicate  self-extinguishing  behavior  in  less  than  30  seconds,  without 
and  with  flaming  droplets,  respectively.  A  rating  of  HB  means  that  the  specimen  burns  in  a 
horizontal  orientation  at  a  rate  which  is  less  than  that  specified  for  a  particular  sample  thickness. 
The  LOI  and  UL  flammability  test  results  are  indicative  of  the  propensity  for  flame  spread  on  a  ton 
polymer  sample  in  the  absence  of  an  external  heat  flux,  but  they  are  poor  predictors  of  full-scale  fire 
performance  in  a  radiant  environment  [16]. 
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The  data  in  table  2  are  ranked  in  descending  order  of  peak  decomposition  temperature  (short-term 
thermal  stability).  High  thermal  stability  results  from  strong  primary  mid  secondary  bonds  and  low 
hydrogen  content  all  of  which  favor  recombination/crosslinking  (charring)  reactions  during  thermal 
degradation  rather  than  the  hydrogen  transfer/termination  reactions  which  lead  to  mass  loss  [98]. 
Excluding  the  halogenated  polymers,  it  is  seen  that  high  thermal  stability  roughly  correlates  with 
flame  resistance. 

In  concluding  this  section  we  note  that  the  thermal  decomposition  temperature  is  a  kinetic  parameter 
which  is  heating  rate  dependent  (c.f.,  equations  53  and  60),  unlike  the  glass  transition  and  melting 
temperatures  which  are  equilibrium  (thermodynamic)  properties.  Thus,  the  decomposition 
temperature  calculated  from  molar  group  contributions  obtained  empirically  from  arbitrarily  slow 
heating  rates  under-predicts  the  polymer  burning  surface  temperature  in  a  fire  where  heating  rates 
may  be  orders  of  magnitude  higher  (c.f.,  figure  1 1  and  equation  17). 

4.  THERMODYNAMICS. 

This  section  addresses  the  intensive,  equilibrium  (rate  independent)  quantities  relevant  to  fire 
behavior  which  are  dependent  on,  and  calculable  from,  the  polymer  chemical  composition  and 
structure.  Included  in  this  category  are  the  char  yield,  heat  of  gasification,  and  chemical  heat  of 
combustion. 

4.1  CHAR  YIELD. 

Char  is  the  carbonaceous  solid  which  remains  after  flaming  combustion  of  the  polymer.  The  char 
yield  is  the  mass  fraction  of  char  based  on  the  original  weight  of  material.  Charring  reduces  the 
volatile  fuel  content  in  a  fire  and  is  thought  to  act  as  a  heat  and  mass  transfer  barrier  which  lowers 
the  chemical  heat  release  rate  in  bench-scale  flaming  combustion  tests  (e.g.,  cone/fire  calorimeters). 
Figure  6  demonstrates  that  the  char  yield  in  a  fire  is  equal  to  the  anaerobic  pyrolysis  residue  at  the 
flaming  surface  temperature.  Thus,  char  formation  takes  place  in  an  oxygen-free  environment 
where  solid-state  oxidation  reactions  are  slow  compared  to  polymer  dissociation  and  gas/char 
formation.  The  equivalence  between  the  char  yield  and  pyrolysis  residue  of  a  material  permits  a 
molecular  interpretation  of  this  important  material  fire  parameter  using  the  large  volume  of 
published  thermogravimetric  data  and  its  correlation  with  chemical  structure. 

Equation  47  shows  that  the  kinetic  fuel  generation  rate  of  a  polymer  during  transient  heating 
decreases  linearly  with  increasing  pyrolysis  residue  or  char  yield.  This  trend  of  decreasing 
flammability  with  increasing  pyrolysis  residue  has  long  been  known  for  small  flame  tests  [99]  and 
the  limiting  oxygen  index  [92, 100]  and  for  this  reason  much  work  has  been  focused  on  empirical 
correlations  between  charring  and  polymer  structure  [91,  92].  Figure  15  shows  the  data  in  table  2 
plotted  to  illustrate  the  well-known  relationship  between  limiting  oxygen  index  (LOI)  and  char  yield 
(p)  for  polymers.  Halogenated  polymers  are  clearly  outliers  from  the  correlation  line  plotted  in 
figure  15,  confirming  a  gas  phase  component  to  the  LOI  measurement. 

Figure  15  shows  that  char  yield  alone  is  not  a  good  predictor  of  the  oxygen  concentration  at  flame 
extinguishment  (LOI).  Moreover,  the  extinguishment  condition  is  a  function  of  both  char  yield  and 
gas  phase  chemistry  as  indicated  by  the  pronounced  effect  of  halogens  on  LOI.  Since  die  flame 
temperature  of  diffusion  flames  increases  with  the  oxygen  concentration,  the  flame  heat  flux  also 
increases.  In  fact,  the  flame  heat  flux  of  burning  polypropylene  increases  in  proportion  to  the 
oxygen  fraction  of  the  environment  [101].  Thus,  a  critical  oxygen  concentration  at  extinguishment 
(LOI)  means  that  a  critical  flame  heat  flux  is  required  to  sustain  steady  burning  of  the  specimen. 
Since  the  UL  94  flammability  test  is  conducted  at  a  fixed  (ambient)  oxygen  fraction,  the  flame  heat 
flux  is  not  an  adjustable  parameter  so  material  properties  alone  will  determine  whether  the  material 
will  sustain  flaming  combustion  after  removal  of  the  ignition  source. 
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Char  Yield,  ^  (%) 


FIGURE  15.  RELATIONSHIP  BETWEEN  LIMITING  OXYGEN  INDEX  AND  CHAR  YIELD 
FOR  HALOGEN  (O)  AND  NONHALOGEN  (•)  POLYMERS  FROM  TABLE  2 

Pyrolysis/char  residue  has  the  character  of  a  thermodynamic  quantity  since  it  depends  only  on 
temperature  and  the  composition  of  the  material  through  the  enthalpy  barriers  to  gas  and  char 
formation,  E  ,  Ec,  in  equation  36.  More  precisely,  char  yield  is  a  statistical  thermodynamic  concept 
wherein  the  total  free  energy  of  the  char  system  at  a  particular  (reference)  temperature  is  the  sum  of 
the  individual  group  contributions  [91,92].  Figure  16  shows  the  equilibrium  residual  mass  fraction 
versus  temperature  for  activation  energy  barriers  AEg_c  =  10,  20,  30,  and  40  kJ/mol  according  to 
equation  36  for  Ag/A,.  =  20. 

Van  Krevelen  [91, 92]  has  devised  a  method  for  calculating  the  pyrolysis  residue  (~  char  yield)  of  a 
polymer  from  its  chemical  composition  and  the  observation  that  the  char  forming  tendency  of 
different  groups  is  additive  and  roughly  proportional  to  the  aromatic  (i.e.,  nonhydrogen)  character 
of  the  group.  Analogous  to  equation  62  for  computing  Tp  from  group  contributions,  the  char  yield 
is  calculated  by  summing  the  char  forming  tendency  per  mole  of  carbon  of  the  chemical  groups, 
Cpr,  i  and  dividing  by  the  molecular  weight  of  the  repeat  unit 


Yc  =  ^«M 


carbon 


i^l  ni  CpT.i 

njMj 


.2  njCpri 

x  100  =  ^ -  x  1200 


The  Cft4  are  the  amount  of  char  per  structural  unit  measured  at  850°C  divided  by  12  (the  atomic 
weight  of  carbon),  i.e.,  the  statistical  amount  of  carbon  equivalents  in  the  char  per  structural  unit  of 
polymer.  Negative  corrections  are  made  for  aliphatic  groups  containing  hydrogen  atoms  in 
proximity  to  char-forming  groups  because  of  the  possibility  for  disproportionation  and  subsequent 
volatilization  of  chain-terminated  fragments  which  are  no  longer  capable  of  crosslinking.  The 
method  is  empirical  and  relatively  simple  to  use  as  illustrated  in  die  following  example  to  predict  the 
char  yield  of  the  fire  resistant,  semicrystalline  thermoplastic  poly(etheretherketone)  (PEEK)  having 
the  chemical  structure 
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The  group  contributions  to  char  formation  for  the  chemical  groups  in  PEEK  are  listed  in  table  3. 


0  200  400  600  800  1000  1200 

Temperature,  T  (°C) 

FIGURE  16.  EQUILIBRIUM  MASS  FRACTION  VERSUS  TEMPERATURE  FOR  ENERGY 
BARRIERS  TO  GASIFICATION  (AEg_c  =  10, 20, 30,  and  40  kJ/mol  according  to  equation  36.) 


TABLE  3.  GROUP  CONTRIBUTIONS  TO  THE  CHAR  RESIDUE  OF 
POLY  (ETHERETHERKETONE)  FROM  REFERENCE  91 


Group 

Quantity 

Mi 

(g/mol) 

CFT,i 

(C-equivalents) 

-<0>- 

3 

76.10 

4 

—o— 

2 

16.00 

0 

O 

— c  — 

1 

28.01 

..  0 

The  char  yield  calculated  with  equation  63  for  poly(etheretherketone) 


Yc(850°C)  = 


_  (3-4)  +  (2-0)  +  0 _ _ 

(3-76.10  g/mol)  +  (2-16.00  g/mol)  +  28.01  g/mol 


12  g/mol  x  100  =  50% 


29 


is  in  good  agreement  with  the  measured  pyrolysis  residue  (char  yield)  of  poly(etheretherketone)  in 
table  2  and  in  the  literature  [102,  103]  at  850°C.  The  char  yield  of  polymers  under  anaerobic 
conditions  is  thus  well  described  using  the  additive  molar  contributions  of  the  individual  groups 
comprising  the  polymer. 

42  HEAT  OF  GASIFICATION. 

In  section  3.2  the  pyrolysis  activation  energy  Ea  was  identified  with  the  molar  enthalpy  of 
gasification  AH*.  Since  E  is  a  thermodynamic  state  function,  it  is  the  sum  of  the  enthalpies  to 
bring  the  polymer  solid  from  the  solid  state  at  the  initial  (room)  temperature  Tc  to  the  gaseous  state 
at  the  decomposition  temperature  Tp.  If  the  stored  heat  is  AHS,  the  enthalpy  of  fusion  (melting)  for 
semicrystalline  polymers  is  AHf,  the  bond  dissociation  enthalpy  is  AH<j,  and  the  enthalpy  of 
vaporization  of  die  decomposition  products  is  AHV ,  then 

Ea  =  AH*  =  AHS  +  AHf  +  AHd  +  AHV  (64) 

Table  4  illustrates  the  magnitude  of  these  enthalpic  terms  on  a  mass  basis  for  amorphous 
poly(methylmethacrylate),  polystyrene,  and  semicrystalline  polyethylene.  Values  are  as  reported  in 
Joules  per  gram  (J/g)  or  have  been  converted  to  a  mass  basis  by  dividing  the  molar  heat  by  the 
molecular  weight  of  the  gaseous  decomposition  products  M„.  The  stored  heat  Ahs  was  obtained  by 
numerical  integration  of  heat  capacity  versus  temperature  [104]  from  ambient  up  to  the  dissociation 
temperature.  The  dissociation  (bond  breaking)  enthalpy  Ahd  is  assumed  to  be  equal  to  the  heat  of 
polymerization  but  opposite  in  sign  for  these  polymers  which  thermally  degrade  by  random  or  end- 
chain  scission  [105].  The  degradation  product  for  polyethylene  is  assumed  to  be  a  tetramer  (i.e., 
octane  with  M  =112  g/mol)  for  the  purpose  of  calculating  the  heats  of  dissociation  and 
vaporization  ona  mass  basis  for  this  polymer,  and  the  degree  of  polyethylene  crystallinity  is  taken 
to  be  90  percent.  All  other  enthalpies  in  table  4  were  obtained  from  handbooks  [104,  106]  using 
individual  monomer  molecular  weights  M  to  convert  the  energies  to  a  mass  basis.  Three  values  of 
h  are  listed  for  each  polymer  in  table  4  for  comparison.  The  uppermost  values  were  calculated  by 
summing  the  enthalpies  listed  in  the  table  for  each  polymer.  The  second  row  of  values  of  hg  are 
calculated  as  Ea/M  from  handbook  values  for  Ea  [104]  and  the  Mg’s  listed  in  the  table.  The  hg 
values  in  the  last  row  of  table  4  are  the  slope  of  mass  loss  rate  versus  external  heat  flux 
measurements  made  in  fire  calorimeters  [46]  according  to  equation  15. 

TABLE  4.  ENTHALPIES  OF  GASIFICATION  FOR  PMMA,  PS,  AND  PE 


Polymer 

PMMA 

PS 

PE 

M  (g/mol) 

100 

104 

28 

M„  (g/mol) 

100 

104 

112 

- &  - 

Value,  J/g 

Ahs 

740 

813 

803 

Ahf 

amorphous 

amorphous 

243 

Ahd 

550 

644 

910 

Ahv 

375 

387 

345 

ha  =  2Ahi 

1665 

1850 

2301 

h„  =  EJM. 

1550  ±  300 

1930 

2230  ±  630 

h„  (measured) 

1600 

1800 

2300 

30 


The  agreement  between  the  enthalpies  of  gasification  calculated  from  the  component  enthalpies 
SAhj  for  the  solid-»gas  phase  change,  the  reported  activation  energy  for  pyrolysis  [107]  Ea/Mg,  and 
the  reported  heat  of  gasification  hg  for  the  three  polymers  in  table  4  is  better  than  ±  10%  when  the 
molar  basis  for  the  activation  energy  is  the  molecular  weight  of  the  degradation  products.  If  Mg  is 
the  average  molecular  weight  of  the  decomposition  products,  then  [44] 


-Ml-  la 

g  "  Mg  "  Mg 


(65) 


and  the  molecular  weights  of  the  decomposition  products  and  the  starting  monomer  should  be  in 
the  ratio 


Mg  Ea 

M  Mhg  (66) 

Polymers  which  pyrolyze  to  monomer  by  end-chain  scission  (depolymerize/unzip)  at  near- 
quantitative  yield  such  as  polymethylmethacrylate,  polyoxymethylene,  and  polystyrene  should  have 
Mg  equal  to  the  monomer  molecular  weight,  M,  i.e.,  Mg/M  ~  1.  Polymers  such  as  polyethylene  and 
polypropylene  which  decompose  by  random  scission  to  multimonomer  fragments  should  have 
Mg/M  >  1.  In  contrast,  polymers  with  high  molecular  weight  repeat  units  (M  >  200  g/mol)  such  as 
nylon,  cellulose,  or  polycarbonate  which  degrade  by  random  scission,  cyclization,  small-molecule 
splitting,  or  chain  stripping  of  pendant  groups  (e.g.,  polyvinylchloride)  will  yield  primarily  low 
molecular  weight  species  (water,  carbon  dioxide,  alkanes,  HC1)  relative  to  the  starting  monomer  and 
should  have  Mg/M  <  1.  Table  5  shows  the  molecular  weight  ratio  Mg/M  calculated  as  Ea/Mhg 
according  to  equation  66  for  some  of  the  commercial  polymers  listed  in  table  2.  Global  pyrolysis 
activation  energies  for  the  thermally  stable  engineering  plastics  listed  in  the  last  four  rows  of  table  5 
are  estimated  to  be  in  the  range  Ea  =  275  ±  25  kJ/mol.  Qualitative  agreement  is  observed  between 
the  modes  of  pyrolysis  (end-chain  scission,  random  scission,  chain  stripping)  and  the  calculated 
fragment  molecular  weight  using  equation  66,  suggesting  that  the  global  pyrolysis  activation  energy 
determined  from  mass  loss  rate  experiments  is  the  molar  enthalpy  of  pyrolysis  of  the  degradation 
products.  The  heat  of  gasification  per  unit  mass  of  original  solid  hg  =  (l-p)Lg  increases  from 
about  2.5  to  3.5  kJ/g  as  a  result  of  the  increase  in  stored  heat  at  the  higher  decomposition 
temperature  of  thermally  stable  polymers  such  as  PC,  PEI,  PPS,  PEEK,  and  PAL 

In  practice,  the  heat  of  gasification  per  unit  mass  of  solid  hg  is  rarely  calculated  because  detailed  and 
reliable  thermodynamic  data  for  the  polymer  and  its  decomposition  products  are  generally 
unavailable  except  for  the  most  common  polymers.  Laboratory  thermogravimetric  analyses  to 
determine  hg  from  the  pyrolysis  activation  energy  E?  require  a  separate  determination  of  the  average 
molecular  weight  of  the  gaseous  thermal  decomposition  products  (equation  65).  The  data  in  table  5 
suggest  that  Mg  ~  75-200  g/mol  for  most  polymers  which  is  in  agreement  with  direct  measurements 
Mg  =  125  ±  50  g/mol  for  the  major  pyrolysis  products  of  aromatic  polyimides  using  gas 
chromatography-mass  spectrometry  [108].  Direct  laboratory  measurement  of  hg  using  differential 
thermal  analysis  (DTA)  [93]  and  differential  scanning  calorimetry  (DSC)  [46]  have  been  reported, 
but  hg  is  usually  measured  indirectly  as  a  component  of  Lg  in  a  constant  heat  flux  gasification 
device  or  fire  calorimeter  [46]. 

4.3  HEAT  OF  COMBUSTION. 

As  stated  previously  polymers  do  not  bum-their  decomposition  products  do.  Thus,  it  is  the  net 
heat  of  complete  combustion  of  the  fuel  gases,  h°c ,  in  equation  16  which  must  be  known  or 
calculable.  It  was  shown  that  polymers  thermally  decompose  to  volatile  fuel  species  and,  possibly,  a 
solid  charred  residue  if  the  molecular  density  of  aromatic/hydrogen-deficient  groups  is  sufficient  to 
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TABLE  5.  HEATS  OF  GASIFICATION,  PYROLYSIS  ACTIVATION  ENERGY,  CHAR 

yield,  and  calculated  molecular  weight  of  decomposition  products 

FOR  SOME  POLYMERS  LISTED  IN  TABLE  2 


Polymer 

M 

(g/mol) 

Lg 

(kJ/g) 

V 

(g/g) 

hg 

(kJ/g) 

Ea 

(kJ/mol) 

Ea/Mhg 

(Mg/M) 

Pyrolysis 

Products 

PP 

42 

2.5 

0 

2.5 

243 

2.3 

C2-C90  saturated  and 
unsaturated 
hydrocarbons 

PE 

28 

2.4 

0 

2.4 

264 

3.9 

PS 

104 

2.2 

0 

2.2 

230 

1.0 

40-60%  monomer 

PMMA 

100 

1.6 

0 

1.6 

160 

1.0 

100%  monomer 

POM 

30 

2.7 

0 

2.7 

84 

1.0 

100%  monomer 

PA  66 

226 

2.6 

0 

2.6 

160 

0.3 

H20  ,  C02 ,  C5  HC  s 

PVC 

62 

2.5 

0.1 

2.8 

110 

0.7 

HC1,  benzene,  toluene 

Cellulose 

162 

3.2 

0.2 

2.6 

200 

0.5 

H20  ,  C02 ,  CO 

PT 

131 

5.0 

0.6 

2.0 

178 

0.3 

Complex  mixture  of 
low  molecular  weight 
products 

PC 

254 

5.3 

0.3 

3.7 

200 

<1 

PEI 

592 

7.7 

0.5 

3.9 

=  275 

<1 

PPS 

108 

6.6 

0.5 

3.3 

=  275 

<1 

PEEK 

288 

5.6 

0.5 

2.8 

=  275 

<1 

PAI 

356 

8.0 

0.6 

3.2 

=  275 

<1 

PX 

180 

6.4 

0.7 

2.2 

=  275 

<1 

allow  inter-  or  intra-molecular  crosslinking  and  cyclization.  The  net  heat  of  complete  combustion 
of  the  polymer  h“  is  related  to  the  net  heats  of  complete  combustion  of  the  volatiles  h°c  and  char 


h  c,n  as 


Kp-lilC 

(1-F) 


(67) 


Thus  the  char  fraction  p  and  its  heat  of  combustion  hCJl  must  be  known  for  charring  materials  to 
calculate  the  heat  of  combustion  of  the  fuel  gases  from  the  heat  of  combustion  of  the  polymer. 
Elemental  analysis  of  the  charred  fire  residue  from  burned  engineering  polymers  [109]  gives  an 
approximate  carbon/hydrogen  molar  ratio  of  5/2  so  that  the  char  repeat  unit  is  C5H2  and 

thermochemistry  calculations  (see  following  section)  give  hc  (1  =  37  kJ/g  which  compares  with  the 
average  polymer  heat  of  combustion  in  table  6,  h°c  p  =  29  ±  6  kJ/g.  For  a  char  yield  of  50% 


(p  =  0.5)  equation  67  predicts 


h! 


»  lh° 

4  “c,p 


(68) 


showing  that  the  heat  of  combustion  of  the  fuel  gases  for  a  charring  polymer  in  a  fire  can  be 
significantly  less  than  the  heat  of  combustion  of  the  polymer  itself. 

At  constant  pressure  and  when  no  nonmechanical  work  is  done,  the  heat  (Q,  q)  and  enthalpy  (H,  h) 
of  a  process  are  equal.  The  flaming  combustion  of  polymers  at  atmospheric  pressure  satisfies  these 
conditions.  The  high-pressure  adiabatic  combustion  of  a  polymer  in  a  bomb  calorimeter  satisfies 
these  conditions  approximately,  since  the  fractional  pressure  change  is  small.  Consequently,  the 
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TABLE  6.  MEASURED  AND  CALCULATED  GROSS  HEATS  (Qc)  OF  POLYMER 
COMBUSTION  (Use  equation  72  to  convert  to  net  heat  of  combustion,  h°  p .  All  values  are  kJ/g.) 


Polymer 

Measured 

From  Oxygen 
Consumption 

From  Heats  of 
Formation 

Polyoxymethylene 

17.39 

15.44 

18.20 

Polytetrafluoroethylene 

8.38 

7.57 

Polyvinylalcohol 

khi 

25.81 

Polyethylene 

47.74 

46.00 

19.53 

18.78 

N/A 

Polypropylene 

46.00 

Polymethylmethacrylate 

26.81 

26.91 

27.50 

Poly(  1 ,4-phenylenesulfide) 

27.98 

Poly(2, 6-dimethyl- 1 ,4  phenyleneoxide) 

34.65 

■ESUSflHif 

Polystyrene 

42.00 

mmmm 

Polyethyleneterephthalate 

heekqhhi 

Epoxy  novolac 

31.37 

36.58 

32.06 

Poly  ( 1 ,4-pheny  leneethersulfone) 

25.42 

23.34 

25.70 

27.91 

25.97 

Poly(hexamethyleneadipamide) 

30.90 

32.74 

Poly(etherketone) 

1 EH 

31.45 

Poly(benzoyl- 1 ,4-phenylene) 

mm 

34.74 

35.90 

Poly(p-phenylene  benzobisoxazole) 

29.18 

Poly(m-phenyleneisophthalamide) 

26.45 

Aramid-arylester  copolymer 

25.27 

28.22 

Poly(p-phenyleneterephthalamide) 

26.92 

28.22 

ksbi i 

Poly(amideimide) 

24.97 

25.28 

26.75 

Poly(acrylonitrile-butadiene-styrene) 

39.43 

Bisphenol  E  Cyanate  Ester 

29.58 

29.62 

Polycarbonate  of  bisphenol-A 

Kim 

30.92 

31.20 

Hexafluorobisphenol  A  Cyanate  Ester 

—  1 1"  m'jL'  ■■■ 

20.00 

19.79 

Bisphenol  A  Cyanate  Ester 

HEEQHMI 

Poly(etheretherketone) 

30.75 

Tetramethylbisphenol  F  Cyanate  Ester 

32.12 

31.72 

Poly(etherketoneketone) 

31.15 

31.50 

22.22 

33.40 

24.72 

29.63 

28.84 

28.83 

Bisphenol  M  Cyanate  Ester 

34.39 

34.15 

33.81 

Polysulfone 

31.20 

Poly(bisphenol-A/aniline)  benzoxazine 

35.45 

35.80 

Polyetherimide 

29.21 

30.00 

Polyester  of  hydroxybenzoic/ 
napthoic  acids 

26.54 

26.22 

26.81 

terms  heat  and  enthalpy  are  used  interchangeably  in  the  following  section  on  methods  for 
determining  the  heat  of  combustion  of  polymers.  Table  6  lists  experimental  values  for  the  gross 
heat  of  combustion  measured  by  oxygen  bomb  calorimetiy  [1 10]  according  to  standard  methods 
[111].  The  second  and  third  columns  in  table  6  list  the  gross  heats  of  combustion  calculated  from 
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oxygen  consumption  thermochemistry  and  from  molar  group  additivity  of  the  heats  of  formation, 
respectively,  according  to  the  methods  described  in  the  following  sections.  Analysis  of  the  data  in 
table  6  shows  that  thermochemical  calculations  based  on  oxygen  consumption  thermochemistry  or 
additive  heats  of  formation  are  within  about  5%  of  measured  values  for  the  gross  heat  of 
combustion  of  polymers. 

4.3.1  Calculation  From  Oxygen  Consumption. 

Heats  of  combustion  of  organic  macro-molecules  can  be  calculated  from  the  oxygen  consumed  in 
the  combustion  reaction  [112].  Oxygen  consumption  is,  in  fact,  the  basis  for  most  modem  bench- 
and  full-scale  measurements  of  heat  release  in  fires  [113,  114].  The  principle  of  oxygen 
consumption  derives  from  the  observation  that  for  a  wide  range  of  organic  compounds,  including 
polymers,  the  heat  of  complete  combustion  per  mole  of  diatomic  oxygen  consumed  is  a  constant,  E, 
which  is  independent  of  the  composition  of  the  polymer.  Mathematically, 


E  = 


nM 
no2  M02 


h° 

=  13.1  ±  0.7  kJ/g-02 


(69) 


where  h°c>p  is  the  net  heat  of  complete  combustion  of  the  polymer  solid  with  all  products  in  their 
gaseous  state,  n  and  M  are  the  number  of  moles  and  molecular  weight  of  the  molecule  or  polymer 
repeat  unit,  respectively,  no2  is  the  number  of  moles  of  O2  consumed  in  the  balanced 

thermochemical  equation,  and  Mo2  =  32  g/mol  is  the  molecular  weight  of  diatomic  oxygen.  In 
equation  69  the  quantity  r0  =  [n02Mo2/nM]  is  the  stoichiometric  oxygen-to-fuel  mass  ratio. 

To  illustrate  the  thermochemical  calculation  of  the  net  heat  of  combustion  we  use  as  an  example 
poly(methylmethacrylate)  (PMMA)  which  has  the  chemical  structure 

CH, 

I  3 

"tCH2-C 

C =0  Poly(methylmethacrylate) 

O 

1 

ch3 

The  methylmethacrylate  repeat  unit  shown  in  brackets  has  the  atomic  composition  C5H802  so  the 
balanced  chemical  equation  for  complete  combustion  is 

C5H802  +  6  02  — >  5  C02  +  4  H20  (70) 

From  equation  70  it  is  seen  that  6  moles  of  02  are  required  to  completely  convert  one  mole  of 
PMMA  repeat  unit  to  carbon  dioxide  and  water.  Inverting  equation  69 


*C,  p 


nQ2  Mq2 


n  M 


(13.1  kJ/g-02)(6  mol  02)(32  g  02/mol  02) 
(1  mol  PMMA)(100  g/mol  PMMA) 


=  25.15  kJ/g  (71) 


To  compare  the  net  heat  of  combustion  to  the  gross  heat  of  combustion  Qc  determined  by  oxygen 
bomb  calorimetry,  the  latent  heat  of  vaporization  of  the  water  produced  in  the  combustion 
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reaction  must  be  added.  An  approximation  based  on  the  weight  fraction  of  hydrogen  in  the 
molecule  wh  is  [115] 


Qc  =  h“p  +  (21.96  kJ/g)wH 


(72) 


For  PMMA,  wh  =  (8  g  hydrogen)/(100  g  monomer)  =  0.08  so  from  equation  72,  Qc  =  25.15  + 
(21.96)(0.08)  =  26.91  kJ/g.  This  is  in  good  agreement  with  literature  values  Qc  =  26.20  kJ/g  [1 16] 
and  26.64  kJ/g  [1 15]  for  PMMA. 


4.3.2  Calculation  From  Heats  of  Formation. 


Calculation  of  the  heat  of  the  combustion  reaction  of  polymers  can  be  carried  out  using  the 
principle  of  molar  additivity  of  the  heats  of  formation  of  the  combustion  products  and  reactants 
[83, 117]  analogous  to  the  previous  calculations  of  decomposition  temperature  and  char  yield  from 
molar  quantities.  The  concept  derives  from  the  fact  that  enthalpy  is  a  state  function  and  therefore  its 
change  in  any  process  is  independent  of  the  path  from  reactants  to  products.  Thus,  the  overall 
enthalpy  of  a  reaction  is  simply  the  sum  of  the  enthalpies  of  the  component  reactions.  In  practice, 
the  heat  of  combustion  of  the  reaction  can  be  calculated  by  subtracting  the  heat  of  formation  of  the 
products  from  the  heat  of  formation  of  the  reactants 


AH”P 


"  np,i  Hp,i  nr,i 


where  p  and  r  denote  products  and  reactants,  respectively,  in  the  standard  state  at  298  K. 


(73) 


For  polymeric  reactants  the  molar  heat  of  formation  can  be  estimated  from  the  tabulated  molar 
contributions  of  the  chemical  groups  which  constitute  the  monomer  or  repeat  unit.  Using  PMMA 
as  an  example  again  with  the  monomer/repeat  unit  chemical  structure  above,  the  heats  of  formation 
of  the  methylmethacrylate  constituent  groups  at  T  =  298  K  are  listed  in  table  7  after  Van  Krevelen 
[91]. 


TABLE  7.  GROUP  CONTRIBUTIONS  TO  THE  HEAT  OF  FORMATION  OF  PMMA 


Group 

Qty. 

Hr>j 

(J/mol) 

nr,|H°i 

(kJ/mol) 

— CH, 

2 

-46,000  +  95T 

-  35.38 

- 1 

-c- 

1 

1 

20,000  +  140T 

61.72 

-CH,- 

1 

22,000 +102T 

8.40 

—O— 

1 

-120,000  +  70T 

-99.14 

0 

II 

— c— 

1 

-132,000  +  40T 

- 120.08 

Total 

-184.48 

Summing  these  group  contributions  gives  the  molar  heat  of  formation  of  the  methylmethacrylate 
monomer,  AHf  =  —184.48  kJ/mol  at  standard  conditions  (T  =  298  K).  The  stoichiometry  of 
complete  combustion  is 

Reactants  — » Products  =  C5H8O2  +  6  O2  — >  5  CO2  +  4  H2O 
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The  tabulated  standard  heats  of  formation  of  the  products  are 

Hp  (H20)  =  -  241.8  kJ/mol;  Hp  (C02)  =  -  393.5  kJ/mol 
The  standard  heats  of  formation  of  the  reactants  are 

Hr  (02)  =  0  kJ/mol;  H°r  (PMMA)  =  -184.5  kJ/mol 

The  molar  heat  of  combustion  of  PMMA  is  therefore 


AHc(PMMA)  = 


Hprod  —  Hreact 

[5  C02  +  4  H20]  -  [C5H802  +  6  02] 

[5  (-  393.5  kJ/mol)  +  4  (-  241.8  kJ/mol)]  -  [-  184.5  kJ/mol  +  6  (0)] 
-2748.7  kJ/mol 


The  absolute  value  of  the  gross  heat  of  combustion  per  unit  mass  is  then 


Qc  (PMMA)  =  |AHC/M| 


=  [2748.7  kJ/mol]/[100  g-MMA/mol] 

=  27.5  kJ/g 

which  compares  favorably  to  literature  values  Qc  =  26.20  kJ/g  [116]  and  26.64  kJ/g  [115]  for 
PMMA. 

5.  THERMAL-CHEMICAL  COUPLING:  THE  PYROLYSIS  ZONE. 

In  this  section  the  solid-state  thermochemistry  of  gaseous  fuel  generation  is  related  to  steady 
macroscopic  burning  through  a  characteristic  dimension-the  pyrolysis  zone  thickness.  It  is  seen 
that  the  pyrolysis  (process)  zone  provides  the  coupling  between  thermal  diffusion  and  solid-state 
thermochemistry.  According  to  equation  52  the  peak  (time-independent)  kinetic  heat  release  rate 
for  complete  combustion  of  volatiles  generated  during  transient  anaerobic  pyrolysis  at  heating  rate 
(3  is 


0c(W/kg) 


lO  fflmax 

nc  m0 


oP(l~l-0Ea 
c  eRTp 


(74) 


From  equation  16  the  steady  macroscopic  heat  release  rate  per  unit  area  of  burning  surface,  S,  for  a 
pyrolysis  zone  of  constant  thickness,  8,  (see  figure  4)  is 


qc(W/m2)  =  Xh°cf  =  XP8h 


O  til 
c  m0 


(75) 


where  mo/S  =  p8  is  the  areal  density  of  pyrolyzing  polymer.  From  equations  74  and  75  the  ratio  of 
the  macroscopic  and  kinetic  heat  release  rates  at  comparable  surface  heating  rates  is 


(76) 
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Thus,  the  macroscopic  and  kinetic  heat  release  rates  are  related  through  the  gas  phase  combustion 
efficiency  %,  polymer  density  p,  and  pyrolysis  zone  thickness  8.  Since  X  and  p  do  not  vary  greatly 
between  materials,  equation  76  shows  that  the  pyrolysis  zone  depth  provides  the  primary  coupling 
between  the  chemical  kinetics  and  thermal  diffusion  in  steady  burning.  Thus,  the  concept  of  a 
distinct  process  zone  is  integral  to  the  solid-state  thermochemistry  of  flaming  combustion. 

The  pyrolysis  zone  thickness  can  be  estimated  using  the  criteria  that  the  mass  loss  rate  falls  to  1/e 
of  the  surface  (x  =  0)  value  at  1/e  of  the  pyrolysis  zone  thickness  x  =  8/e.  Since  the  mass  loss  rate 
has  an  Arrhenius  form 


m(8/e)  1 

Ea 

1  iV 

— -  =  —  =  exp 

m(0)  e 

R 

[T(8/e)  Tp  J 

(77) 


Equation  77  shows  that  the  mass  loss  rate  to  falls  to  1/e  of  the  surface  value  for  typical  Ea  =  200  ± 
50  kJ/mol  when  T(8/e)  is  10-20  K  lower  than  the  surface  temperature  Tp,  i.e.,  Tp  -  T(8/e)  =  20  K. 
Since  Tp  -  T(8/e) «  Tp  «  T(8/e),  equation  77  takes  the  simple  form 


Tp  -  T(8/e) 


(78) 


Substituting  equation  1 1  into  equation  10  with  x  =  8/e 

Tp-T(8fe)  =  (Tp-T0)(l-exp[-&g] 

Eliminating  Tp  -  T(8/e)  between  equations  78  and  79  for  RTp2/(Ea(Tp  -  T0)) «  1, 


(79) 


v  8^ 
a  e 


rt: 


(80) 


or  with  equation  1 1 


s  K  eRTp 

8  =  - - ^ — 

9  net  a 


(81) 


For  typical  polymer  values  Tp  =  750  K,  Ea  =  200  kJ/mol,  and  K(Tp)  «  0.2  W/m-K,  equation  81 
predicts  8  ~  0.3  mm  at  a  net  incident  heat  flux,  qne,  =  50  kw/m2,  which  is  in  agreement  with 
estimates  8  <  1  mm  [32]. 

A  rate-independent  material  flammability  parameter  emerges  from  this  analysis  when  the  kinetic 
heat  release  rate  on  the  right  side  of  equation  74  is  normalized  for  heating  rate 


ric(J/g-K)  ^ 


P 


h°c(l-p)Ea 

eRT„2 


(82) 
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The  thermokinetic  flammability  parameter  1%  has  the  units  and  significance  of  a  heat  [release] 
capacity  (J/g-K)  when  the  linear  heating  rate  is  P(K/s).  Substituting  the  heat  release  capacity  T]c 
(equation  82),  pyrolysis  zone  depth  8  (equation  81),  and  the  heating  rate  (3  (equation  17)  into 
equation  75  recovers  the  thermal  diffusion-limited  steady  heat  release  rate  (equation  16)  from  the 
thermochemistry  after  cancellation  of  terms 


Ac 


XP5Pric  =  XP 


k  eRTp 

•  2  ' 
Qnet 

h°  (1  ~  M)Ea 

Qnet 

Kphg 

eRTp 

(83) 


Obtaining  the  macroscopic  heat  release  rate  from  the  derived  diffusion  (P)  thermokinetic  (r|c)  and 
thermal-chemical  coupling  (8)  parameters  shows  that  the  thermochemistry  of  diffusion-limited 
steady  burning  is  self-consistent. 

6.  CHEMICAL  STRUCTURE  AND  FIRE  BEHAVIOR. 

Sections  1  through  5  demonstrate  that  the  time  to  ignition  and  steady  burning  rate  of  a  polymer  in  a 
particular  fire  environment  (i.e.,  external  heat  flux)  can  be  expressed  in  terms  of  material  properties 
which  are  calculable  from  the  chemical  structure  of  the  polymer  and/or  measurable  using  standard 
laboratory  thermal  analyses.  In  this  concluding  section  we  apply  chemical  structure-fire  property 
relationships  to  the  calculation  of  polymer  flammability,  heat  flux  dependence  of  the  burning 
temperature,  and  heat  release  rate. 

6  1  FLAMMABILITY. 

The  data  in  figure  15  and  table  2  suggest  that  both  the  halogen  content  and  the  char  yield  \x  are 
needed  to  correlate  all  of  the  polymer  flammability  data.  If  the  effect  of  halogens  on  polymer 
flammability  is  primarily  to  alter  the  gas  phase  combustion  efficiency  %  [118]  the  concept  of  a 
critical  heat  release  rate  developed  in  section  2.2  can  be  extended  to  flame  extinguishment  in  the  UL 
94  test. 

Figure  17  is  a  plot  of  UL  94  ranking  versus  the  quantity  |i+LOI  for  all  of  the  polymers  in  table  2. 
It  is  apparent  from  this  plot  that  a  p+LOI  value  below  about  55%  is  associated  with  sustained 
burning  (UL  HB)  while  polymers  with  p+LOI  greater  than  55%  tend  to  self-extinguish 
(UL  V-0/1/2)  after  removal  of  the  small  flame  ignition  source.  Polymers  with  values  of  p+LOI  > 
95%  do  not  bum  or  ignite  even  after  repeated  attempts  at  ignition  (e.g.,  UL  94-5V)  and  also  exhibit 
very  low  heat  release  rate  in  bench-scale  fire  tests  [16,  17,  18].  These  observations  are  consistent 
with  a  critical  energy  release  rate  (section  2.2)  for  burning  in  the  UL  94  test. 

Combining  equations  9  and  16  the  energy  release  rate  in  the  UL  94  test  can  be  represented  as 

h°c  [•  •  •  \ 

qc  -  X  (1  _  lO  A  external  ^  Aflame  A  critical  j 

Accounting  for  gas  phase  combustion  inhibition  with  the  approximation  X  *  (J  -  L01)  derived 
from  combustion  efficiency  data  [46]  and  limiting  oxygen  indices  [56]  of  burning  polymers,  the 
non-burning  condition  after  removal  of  the  Bunsen  burner  flame  in  the  UL  94  test  when  qext  =  0  is 
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H  +  LOI 


1 


qc,c 


=  l- 
«  0.56 


.  v 

(  3  kJ/g  \  t 
^30kJ/gj 


Aflame  ^critical  J 


+  (p-LOI) 


100  kW/m 


35  kW/m2  - 15  kW/m2 


+  (0.25)(0.25) 


where  typical  values  for  the  heat  of  gasification  (table  5),  heat  of  combustion  (table  6),  flame  heat 
flux  [119],  critical  heat  flux  [46],  and  the  critical  energy  release  rate  qc,cr  =  100  kW/m2 
(section  2.2)  were  assumed.  The  small  product  term  p,LOI  was  evaluated  at  |i  =  LOI  =  0.25 
corresponding  to  approximate  values  at  the  flammability  transition  p+LOI  «  0.5  in  figure  17.  The 
critical  energy  release  rate  criteria  for  burning  in  the  UL  94  test  predicts  a  transition  (V-l,  V-2)  from 
flammable  (HB)  to  nonflammable  (V-0)  behavior  for  both  halogen  and  nonhalogen  polymers  at 
p+LOI  =  0.56.  This  prediction  is  in  excellent  agreement  with  the  data  in  figure  17  and  suggests  a 
predominantly  additive  effect  of  solid-state  and  gas-phase  mechanisms  of  flame  resistance  (i.e., 
p+LOI)  with  synergistic  effects  embodied  in  the  product  term  p,-LOI. 


FIGURE  17.  UNDERWRITERS  LABORATORY  UL  94  RANKING  OF  FLAMMABILITY 
VERSUS  THE  COMBINED  PARAMETER  p+LOI  FOR  THE  POLYMERS  IN  TABLE  2 
SHOWING  TRANSITION  FROM  BURNING  TO  SELF-EXTINGUISHING  BEHAVIOR  AT 

p+LOI «  0.56  (56%) 


6.2  BURNING  TEMPERTURE. 

In  section  2  we  made  the  ad  hoc  assumption  that  the  ignition  and  steady-burning  temperatures  at 
constant  heat  flux  in  the  fire  environment  were  equal  to  the  peak  pyrolysis  temperature  calculated 
from  group  contributions  or  measured  at  constant  heating  rate  in  a  laboratory  environment  (e.g., 
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thermal  analysis  instrument)  in  order  to  establish  a  relationship  between  the  chemical  structure  of 
polymers  and  their  fire  behavior.  To  test  this  hypothesis  we  calculate  the  surface  burning 
temperature  of  PMMA  as  a  function  of  external  heat  flux.  Equation  17  gives  the  effective  heating 
rate  at  the  polymer  surface  which  when  substituted  into  equation  60,  gives  the  peak  pyrolysis 
temperature  as  a  function  of  heat  flux.  Figure  18  shows  data  (private  communication  M.A. 
Delichatsios,  Factory  Mutual  Research,  1999)  for  the  measured  surface  temperature  of  1-inch-thick 
PMMA  slabs  during  steady  nonflaming  gasification  (under  nitrogen)  at  the  indicated  radiant 
external  heat  fluxes.  For  nonflaming  gasification,  qflame  =  0,  so  from  equation  17  with  qext  =  qnet 
»  qcr  and  Kpc  =  5  x  105  W-s-m^-K'2  we  calculate  the  surface  heating  rates  P  shown  along  the  top 
of  figure  18.  From  the  calculated  P  with  p*  =  1013  K/s  (see  figure  14)  and  Ea  =  160  U/mol  for 
PMMA,  equation  60  predicts  the  steady  surface  temperature  versus  heat  flux  relationship  indicated 
by  the  solid  line  in  figure  18.  The  measured  and  calculated  surface  temperatures  for  steady-state 
gasification  are  in  reasonable  agreement  for  PMMA— suggesting  that  thermal-diffusion  limited 
chemical  kinetics  influences  the  fuel  generation  process  for  thermally  thick  samples. 

Surface  Heating  Rate  (K/s) 


1.7  3.6  6.3  12.3  15.7  20.6 


FIGURE  18.  SURFACE  TEMPERATURE  OF  PMMA  VERSUS  EXTERNAL  RADIANT 
HEAT  FLUX  FOR  STEADY  NONFLAMING  GASIFICATION 

63  HEAT  RELEASE  RATE. 

Figure  19  shows  experimental  data  [16,  17,  18,  19,  48,  120,  121]  for  the  3-minute  average  heat 
release  rate  measured  in  a  fire  calorimeter  operating  on  the  oxygen  consumption  principle  [122]. 

Incident  heat  flux  on  the  thermally  thick  polymer  samples  was  q  ext  —  50  kW/m2.  The  trend  of 
decreasing  heat  release  rate  in  flaming  combustion  with  increasing  char  yield  of  the  polymer  is 
obvious  in  figure  19.  However,  the  nonlinear  relationship  between  heat  release  rate  and  char  yield 
suggests  that  simple  rule-of-mixtures  behavior  based  on  the  reduced  fuel  fraction  is  not  the 
underlying  mechanism.  As  discussed  in  section  2  char  acts  as  a  thermal-  and  mass-diffusion 
barrier  as  it  accumulates  on  the  surface  of  burning  polymers.  Table  2  shows  that  charring  polymers 
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also  tend  to  be  more  thermally  stable  and  their  flammability  depends  on  both  solid-state  (charring) 
and  gas  phase  (combustion)  reactions.  Deviation  of  polyvinylchloride  (PVC)  from  the  trend  line 
for  hydrocarbon  polymers  in  figure  19  highlights  the  effect  of  gas  phase  combustion  inhibition  by 
halogens  on  heat  release  rate  as  well  as  smale-scale  flammability  tests. 


FIGURE  19.  AVERAGE  FLAMING  HEAT  RELEASE  RATE  VERSUS  CHAR  YIELD  FOR 
POLYMERS  LISTED  IN  TABLE  2  (Not  listed  in  table  2  are  polypropylene  (PP),  (Polyester)-, 
(Vinylester)-,  and  (Epoxy)-thermosets,  polybutyleneterephthalate  (PBT),  poly(2,6, 
dimethyl)phenlyene  oxide  (PPO),  and  polyethersulfone  (PESF).) 

Section  2  provided  the  relationship  between  the  thermal  (k,  p,  and  c)  and  combustion  (Tign,  |X,  hg, 

and  h°cp )  properties  of  materials  and  their  steady  heat  release  rate  during  flaming  combustion 

(equation  16)  derived  from  a  mass  and  energy  balance  at  the  surface 

•  .  „  •  n  h°c  .  h°c  . 

qc=  %hcmg  =  X(l-H)i-qnet  =  %T-q  net 

ng  •  g 

Methods  were  presented  for  obtaining  the  solid  combustion  properties  from  molar  group  additivity 
calculations  (Tign  «  T ,  p,  and,  h°c  p )  or  thermochemical  calculations  (hg).  Methods  for  calculating 
equilibrium  thermophysical  properties  of  polymers  (k,  p,  c)  and  their  temperature  dependence  from 
molar  group  additivity  calculations  are  well  established  [83,  85],  as  are  the  measurement  techniques 
for  these  properties  [93,  123]  and  were  not  reviewed  in  the  chapter.  Thus,  in  principle  all  of  the 
fundamental  thermal  and  combustion  properties  necessary  to  predict  the  flaming  heat  release  rate  of 
burning  polymers  in  a  given  fire  environment  (%,  qnet )  are  calculable  to  a  first  approximation  from 
the  polymer  chemical  structure.  As  an  example  we  calculate  the  steady  (average)  heat  release  rate  of 
a  char  forming  polymer  (bisphenol-A  polycarbonate)  from  its  chemical  structure  using  the  concepts 
developed  in  this  chapter  and  compare  the  result  to  experimental  data  for  this  polymer  plotted  in 
figure  19. 
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The  chemical  structure  of  bisphenol-A  polycarbonate  (e.g.,  LEXAN™)  is 


CH 


BPA- 

Polycarbonate 


From  the  chemical  groups  in  polycarbonate  we  calculate  the  following  properties  from  their  molar 
group  contributions  [91]: 


» Density  at  pyrolysis  temperature  p  (Tp)  ~  1040  kg/m3 
» Pyrolysis  temperature  Tp  =  683  K 

•  Char  yield  p  =  0.24 

» Net  Heat  of  Combustion  (solid)  h°p  =  29.7  kJ/g 


reference  73 
equation  62 
equation  64 

equations  71  and  72 


Note  that  the  calculated  pyrolysis  temperature  of  683  K  is  significantly  below  the  measured 
pyrolysis  temperature  for  polycarbonate  reported  both  in  table  2  and  reference  52.  From  these 
properties  and  the  literature  value  Ea  =  200  kJ/mol  for  polycarbonate  [124]  we  calculate  the 
following  parameters: 

•  The  net  heat  of  complete  combustion  of  the  fuel  gases 

io  _  hc,p~PhCin  _  29.7  kJ/g  -  (0.24)  37  kJ/g  _  27  4  kJ/g  (equation  67) 
c  “  (l-p)  "  (1-0.24) 

•  The  surface  heating  rate  for  polycarbonate  having  Kpc  =  7.5  x  105  W-s  ’-m^-K 2  [46] 


p  -  i 


The  pyrolysis  zone  depth 


7.5  x  10 


(50kW/m2)2 

5_^=lV683K-300k) 

m-K2  1  1 


»  4.4  K/s 


(equation  17b) 


8  =  - 


eK  RT2  (2.718)(02  W/m-K)(8.314  J/mol-K)(683  K) 


(50  kW/m2)(200  kJ/mol) 


(equation  81) 


=  2  x  10-4  m 


The  heat  release  capacity 


nc(j/g-K)  ^ 


h°c(l-p)Ea 


(equation  82,  table  5) 


(27.4  kJ/g)(l  -  0.24)(200  kJ/mol) 
(2.7 1 8)(8.3 14  J/mol-K)(683  K)2 


=  395  J/g-K 
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Collecting  these  terms  in  equation  84  and  approximating  the  gas  phase  combustion  efficiency  as 

X  =  1  -  LOI  =  (1  -  0.26)  =  0.74  (see  section  4.1) 

calculate  the  average  heat  release  rate  in  flaming  combustion  from  equation  83 

qc  =  xpSpri  c  =  (0.74)(1040  kg/m3)(2  x  10-4m)(4.4  K/s)(395  J/kg-K) 

=  268  kW/m2 

which  agrees  with  the  experimental  value  qc  =  250  kW/m2  reported  in  figure  19  for  bisphenol-A 
polycarbonate  to  within  the  accuracy  of  the  calculation  and  the  experimental  accuracy  (±15%)  of 
the  heat  release  rate  measurement.  Note  that  the  much  simpler  calculation  using  equation  16  with 
an  average  enthalpy  of  gasification,  h  «  3.0  ±  0.8  kJ/g  (table  5)  also  gives  reasonable  agreement 
with  measured  values  for  the  chemical  heat  release  rate  of  burning  polycarbonate. 

i  ■  (a74>((3kj^a-o.24))(50kw/m2>  "  257kW/m2 
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